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. Abstract 

We study the point spectrum of the nonlinear Dirac equation linearized at one of the solitary wave solutions 
w (a;)e _ . We prove that, in any dimension, the linearized equation has no embedded eigenvalues in the part of the 
essential spectrum beyond the embedded thresholds (located at A = ±i(m + |w|)). We then prove that the birth of 
point eigenvalues with nonzero real part from the essential spectrum is only possible from the embedded eigenvalues, 
and therefore can not take place beyond the embedded thresholds. We also study the birth of point eigenvalues in the 
nonrelativistic limit, oj — y rn. 

We apply our results to show that small amplitude solitary waves (u < m) of cubic nonlinear Dirac equation in 
■ one dimension are spectrally stable. 

-I— > 1 

1 Introduction 

While a lot is known about the nonlinear Schrodinger and Klein-Gordon equations (see e.g. the review [Str89|), there 
1 are still numerous open questions for systems with Hamiltonians that are not sign-definite, such as the Dirac-Maxwell 

^sO . system [Gro66 Wak66| and the nonlinear Dirac equation [Sol70|. There has been an enormous body of research 

devoted to the nonlinear Dirac equation, which we do not hope to cover comprehensively, only giving a very brief 
sketch. The existence of standing waves in the nonlinear Dirac equation was studied in |Sol70, CV86, Mer88 ES95|. 
Local and global well-posedness of the nonlinear Dirac equation was further addressed in [EV97| (semilinear Dirac 
equation in 3D) and in [MNNO05I (nonlinear Dirac equation in 3D). There are many results on the local and global 
well-posedness in ID; we mention [ST10 MNTlOIICanl lllPell 11 . There were several attempts to study the stability 
of solitary waves of the nonlinear Dirac equation analyzing whether the energy functional is minimized with respect to 
dilations and other families of perturbations ]Bog79| [SV861 ICKMS 1 01 |MQC + 1 2| . The spectrum of the linearization 
at solitary waves of the nonlinear Dirac equation in ID was computed in |Chu07, BC12|, suggesting the absence of 
eigenvalues with positive real part for linearizations at small amplitude solitary waves; we will say that such solitary 
waves are linearly stable, or spectrally stable. The numerical simulations of the evolution of perturbed solitary waves 



|MQC + 12| suggest that the small amplitude solitary waves in ID nonlinear Dirac equation are dynamically stable 



(or nonlinearly stable). The asymptotic stability of small amplitude solitary waves in the external potential has been 
studied in [Bou06 Bou08 PS 10 1. The first approach to the translation-invariant case in 3D (based on the spectral 
stability assumptions) is in MBC1 II . 

In the present paper, we study the spectral stability of solitary waves in the nonlinear Dirac equation. More 
precisely, we study the scenarios of the emergence of positive-real-part eigenvalues in the spectrum of the linearization 
at different solitary waves. We use our results to prove that the small amplitude solitary waves to the cubic nonlinear 
Dirac equation in ID are spectrally stable; this is the first rigorous result of this type. 

We consider the nonlinear Dirac equation in K", n > 1: 

m t '>l> = D n d-fWM)M, ^(x,t)eC N , xew\ (l.i) 

where N is even, / € C°°(R), /(0) = 0, and D m is the free Dirac operator: 

D m = -ia • V + j8m, m > 0, (1.2) 
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a = (a 3 ) i<j<n, where a? and (3 are self-adjoint NxN Dirac matrices. See Section|2]for the details. We are interested 
in the stability properties of the solitary wave solutions to ( 11.11 ): 

4>(x, t) = (j> u {x)e-^\ ^ € H°°(R n , C N ). 

We consider the perturbation of a solitary wave, ((^^(x) + p(x, t))e~ lult , where p is a perturbation, and study the 
spectrum of the linearized equation on p. If the spectrum contains eigenvalues with positive real part, then the solitary 
wave is called linearly unstable, and one expects to be able to prove that this linear instability leads to the orbital 
instability, in the sense of |GSS87|. 

If the spectrum of the linearized equation is on the imaginary axis, we will say that the corresponding solitary wave 
is spectrally stable, or linearly stable. In this case, one hopes to prove the asymptotic stability of solitary waves using 
the dispersive estimates similarly to how this has been done for the nonlinear Schrodinger equation. First results in 
this direction are already appearing [BC1 1 1, with the assumptions on the spectrum of the linearized equation playing a 
crucial role. (Note that in the context of the nonlinear Dirac equation we do not know how to prove the orbital stability 
|GSS87| except via proving the asymptotic stability first.) 

Since the isolated eigenvalues depend continuously on the perturbation, it is convenient to trace the behavior of 
"unstable" eigenvalues (eigenvalues with positive real part) for linearization at the solitary waves w e~ lwt , considering 
w as a parameter. For example, if one knows that solitary waves with w in a certain interval are spectrally stable, one 
wants to know how and when the "unstable" eigenvalues may emerge from the imaginary axis. The emergence of 
unstable eigenvalues from A — is described by the Vakhitov-Kolokolov stability criterion [Comll]. In the present 
paper, we investigate the bifurcation of eigenvalues from the essential spectrum. We prove that such bifurcations are 
only possible from the part of the essential spectrum between the edge of the spectral gap and the embedded threshold: 
A G iR., |A| G [m — m + We will show that such bifurcations, if they exist, can only originate from the 
embedded eigenvalues (with the exception of u G {0; ±m}, A = ±i(m — |w|)). 

Our approach to the spectral stability of solitary waves in the nonlinear Dirac equation is also applicable to the 
Dirac-Maxwell system. Let us mention that the local well-posedness of the Dirac-Maxwell system was proved in 
|Bou96|, while the existence of standing waves in the Dirac-Maxwell system is proved in |EGS96| (for uj e (— m, 0)) 
and BAbe98l (for uj € (—m, m)). (For an overview of these results, see MES02I .) According to the nonrelativistic 
asymptotics (uj > — m) from [CS12|, we expect that the solitary waves with u sufficiently close to — m are spectrally 
stable. 

Another situation where our methods are applicable is the analysis of stability of gap solitons in nonlinear coupled- 
mode equations. Such systems appear in the context of photonic crystals [dSS94], where they describe counter- 
propagating light waves interacting with a linear grating in optical waveguides made of material with periodically 
changing refractive index |dSSS96 GWH01 ]. Coupled-mode systems also describe matter-wave Bose-Einstein con- 
densates trapped in an optical lattice |PSK04|. The numerical analysis of the spectrum of the linearizations at the gap 
solitons is performed in |BPZ98, CP06]. The stability analysis of small-amplitude gap solitons based on the study of 
bifurcations from the embedded eigenvalues of the linear equation in the external potential is in BGW08L 

Derrick's theorem 

As a warm-up, let us consider the linear instability of stationary solutions to a nonlinear wave equation, 

-i> = -Aip + g(il>), ip = ip{x,t) € R, xeR n , n>l, tel. (1.3) 

We assume that the nonlinearity g(rf), rj 6 R is smooth. Equation ( 11.31 ) can be written as a Hamiltonian system 
7T = —S^E, ijj = 5 n E, with the Hamiltonian E{ijj, tt) = J Rn + - + G(V>)) dx, where G(rj) = Q g(Q d£. 
There is a well-known result about non-existence of stable localized stationary solutions, known as Derrick's theorem: 

Lemma 1.1 (Derrick's theorem [Der64|). Equation (|1.31 l can not have stable, time-independent, localized solutions 
in three dimensions. 

Here is the argument from [Der64|, which applies to dimensions n > 3. Denote 

T(9) = - I \\79\ 2 dx, V(d)= I G{6)dx. 
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lfifj(x, t) = 9{x) is a localized stationary solution, so that = ip ■ 
e x {x) = 6{x/X), using the identities T{6 X ) = A"- 2 T(6>), V{6 X ) 



$f(M),0 = 7r = 
X n V(0), one has: 



SE i 
' Sip 1 



, 0), then for the family 



= (^(M), |£k=i) = 5aL =1 ^a,0) = (n - 2)T(0) + nV(0), 



00* 



and it follows that 



dl\ x=1 E(0 x ) = (n - 2)(n - 3)T(0) + n(n - = -2(n - 2)T(0), 



which is negative as long as n > 3. That is, S 2 E < for a variation corresponding to the uniform stretching, and 
the solution 6(x) from the physical point of view is to be unstable. We remark that the fact that d\E(9\)\ x=1 was not 
negative for n — 1 and 2 does not prove that in these dimensions the localized stationary solutions are stable; it just 
means that a particular family of perturbations failed to catch the unstable direction. Let us modify Derrick's argument 
to show the linear instability of stationary solutions in any dimension. 

Lemma 1.2 (Derrick's theorem for n > 1). For any n > 1, a smooth finite energy stationary solution 9 G H°°(M. n ) 
to the nonlinear wave equation (11.3b is linearly unstable. 

Proof. Since 9 satisfies -A6 + g{9) = 0, we also have -Ad Xl 9 + g'(9)d Xl 9 = 0. Due to lim 9(x) = 0, d Xl 9 

\x\— ►oo 

vanishes somewhere. According to the minimum principle, there is a nowhere vanishing smooth function \ <= L 2 (W l ) 
(X € H°°(M n ) due to A being elliptic) which corresponds to some smaller (hence negative) eigenvalue of I = 
—A + g'(9), l\ — —<?X, with c > 0. Taking ip(x, t) — 9(x) + r(x, t), we obtain the linearization at 9, —r = — Ir, 
which we rewrite as 



The matrix in the right-hand side has eigenvectors 



r 




' 1 




r 


s 




-I 




s 



X 



, corresponding to the eigenvalues ±c G R; thus, the solution 



' is linearly unstable. 

Let us also mention that d 2 \ T=0 E(9 + t\) < 0, showing that 5 2 E{9) is not positive-definite. 



□ 



Remark 1.1. A more general result on linear instability of stationary solutions to ( 11.3b is proved in [KS07 |. In par- 
ticular, it is shown there that the linearization at a stationary solution may be spectrally stable when this particular 
stationary solution is not from H 1 (such examples exist in higher dimensions). 

One can see that the linear stability analysis (Lemma [L2l ) is more conclusive than the analysis of whether the energy 
functional is minimized under a particular families of perturbations or not (Lemma fTTl i; in particular, we have just seen 
that, in the context of the nonlinear wave equation (11.3b . the dilation perturbation fails to pick up the unstable directions 
in dimensions n < 2. Moreover, for the nonlinear Dirac equation, it is easy to demonstrate that the solitary waves never 
correspond to the energy minima under the charge constraint, and, although we know that the energy is minimized 
at a solitary wave under some particular charge-preserving perturbations [ |Bog79| TSV86I ICKMSIOI MQC + 12| , it is 
not clear whether any conclusions could be drawn from this. This suggests that in order to have an insight about the 
stability of localized spinor solutions, one needs to perform the linear stability analysis. 



Vakhitov-Kolokolov stability criterion for the nonlinear Schrodinger equation 

For the nonlinear Schrodinger equation and several similar models, real eigenvalues could only emerge from the origin, 
and this emergence is controlled by the Vakhitov-Kolokolov stability condition HVK73I . Let us give the essence of the 
linear stability analysis on the example of the (generalized) nonlinear Schrodinger equation, 

= — — /(l^l 2 )^, ^ = V(M)eC, iel", n>l, teR, (1.4) 

where f(rf) is a smooth function with /(0) = 0. One can easily construct solitary wave solutions (f>(x)e~ lujt , for some 
ui G IR and <f> G H 1 (W l ): <j){x) satisfies the stationary equation ui(f> = —^A<fi — /(0 2 )0, and can be chosen strictly 
positive, even, and monotonically decaying away from x — 0. The value of ui can not exceed 0. We will only consider 
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the case u < 0. We use the Ansatz ip(x, t) - 
is called the linearization at a solitary wave: 



(c/>(x) + p(x, t))e with p(x, t) £ C. The linearized equation on p 

(1-5) 



d tP =\{- ±Ap -up- f{<f> 2 )p - 2/'(0 2 )0 2 Rep 



Remark 1.2. Because of the term with Re p, the operator in the right-hand side is R-linear but not C-linear. 

To study the spectrum of the operator in the right-hand side of ( 11.5b . we first write it in the C-linear form, consid- 

Bep(x,t) 



ering its action onto p(x, t) = 



lmp(x,t) 

ftp = jlp, 



p(x,t) = 



Bep(x 7 t) 
Im p(x, t) 



where 



1 



1 

-1 



I 




l_ 



with l_ 



2 



/(</> 2 ), l+=l_-2/'( 



(1-6) 



(1.7) 



If 4> S <S(R ra ), then by Weyl's theorem on the essential spectrum one has 

Ccss(l-) = cr eS8 (l+) = [|w|,+oo). 

Lemma 1.3. cr(jl) c K U iR. 



Proof. We consider (jl) 2 = — » + . " . Since l_ is positive-definite (<f> G kerl_, being nowhere zero, corre 

(J L_|_ L_ 

sponds to the smallest eigenvalue), we can define the self-adjoint root of l_; then 

^(0l) 2 )\{0} = <7 rf (l_l+)\{0} = a d (l + l_)\{0} = a d (l V2 l + l V2 )\{0} C R, 
with the inclusion due to l]/ 2 \-\-\]/ 2 being self-adjoint. Thus, any eigenvalue A <E <Jd(jV) satisfies A 2 € R. □ 

Given the family of solitary waves, cj) UJ (x)e^ lu ' t , w e O C 1, we would like to know at which cj the eigenvalues 
of the linearized equation with Rc A > appear. Since A 2 6 R, such eigenvalues can only be located on the real axis, 
having bifurcated from A = 0. One can check that A = belongs to the discrete spectrum of jl, with 







= 0, 











for all to which correspond to solitary waves. Thus, if we will restrict our attention to functions which are spherically 
symmetric in x, the dimension of the generalized null space of jl is at least two. Hence, the bifurcation follows the 
jump in the dimension of the generalized null space of jl. Such a jump happens at a particular value of u if one can 





solve the equation jla 



This leads to the condition that 







is orthogonal to the null space of the 



adjoint to jl, which contains the vector 







ft 



2 /2 = 0. A slightly more careful 



; this results in (0 w ,oL0 w / - ^n^n L 2 

analysis [CP03| based on construction of the moving frame in the generalized eigenspace of A = shows that there 
are two real eigenvalues ±A G R that have emerged from A = when lu is such that 9 W 1 1 0^, 1 1 2 , becomes positive, 
leading to a linear instability of the corresponding solitary wave. The opposite condition, 



ft 



<o, 



(1.8) 



is the Vakhitov-Kolokolov stability criterion which guarantees the absence of nonzero real eigenvalues for the nonlin- 
ear Schrodinger equation. It appeared in [ VK73 , CL82 , Sha83 Wei86 GSS87 1 in relation to linear and orbital stability 
of solitary waves. 

The above approach fails for the nonlinear Dirac equation since l_ is no longer positive-definite. Now one no 
longer knows whether the eigenvalues are only on the real or imaginary axes, neither one knows whether ( 11.8b or its 
opposite is needed for stability. All we know is that d 1 -8b vanishes when the eigenvalues collide at A = 0. 

Our conclusions: 
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1. Point eigenvalues of the linearized Dirac equation may bifurcate (as ui changes) from the origin, when the 
dimension of the generalized null space jumps up (at the values of uj when c^H^wlU 2 = 0). 

2. Since the spectrum of the linearization does not have to be a subset of RU iR, there may also be point eigenvalues 
which bifurcate from the imaginary axis (either from the essential spectrum or from the collision of eigenvalues 
in the spectral gap) into the complex plane. We do not know particular examples of such behavior for the 
nonlinear Dirac equation. 

3. Moreover, there may be point eigenvalues already present in the spectra of linearizations at arbitrarily small 
solitary waves. Formally, we could say that these eigenvalues bifurcate from the essential spectrum of the free 
Dirac equation, which can be considered as the linearization of the nonlinear Dirac equation at the zero solitary 
wave. 

The first scenario has been studied in MComl II . Here we will investigate the second and the third scenarios. 

Here is the plan of the paper. The results are stated in Section [2] We consider the properties of solitary waves in 
Section[3] General properties of the linearization at solitary waves are considered in Section|4] Properties of embedded 
eigenvalues are in Section|5] Bifurcations of eigenvalues from the imaginary axis are considered in Sections[6][7] and|8] 
The Hardy-type inequalities and Carleman-Berthier-Georgescu inequalities are considered in AppendiceslAlandlBl 

2 Main results 

Let Q be the operator of multiplication by x and (Q) be the operator of multiplication by y/l + x 2 . 
Foru G L 2 (R"), we denote ||it|| = ||m||l 2 - For s, k € R, we define 

H k s {M, n ) = {ue S'(R n ), \\u\\ Hi < ex)} , ||«|| H j = \\(Q) s (V) k u\\ L 2. (2.1) 

where if k = 0, we write L 2 instead of H®. 
Let 

D m = —ia ■ V + /3m, m > 

be the free Dirac operator. Here a ■ V = X)"=i aJ ~(fb> with a J and /3 being self-adjoint N x N Dirac matrices which 
satisfy 

(a 3 ) 2 = (3 2 = I N , a>a k + a k a 3 = 28 ]k Iff, a 3 (5 + (3a J = 0, l<],k<n. 

Above, Im is the N x N identity matrix. 

We consider the nonlinear Dirac equation ( 11.11 ), 

\Btil> = D m il>-fWPil>)M, iP{x,t)€C N , i£R n , (2.2) 

where / € C°°(R), /(0) = 0. 

Assumption 1. There exists a nonempty open interval O C (— m, m) such that the equation 

uiu = D m u — f(u*j3u)/3u (2.3) 
admits a continuous family of solutions O — > if°°(R n ), u> H> tf> u . 

Then <j) ul (x)e~~ lu ' t is a solitary wave solution to the nonlinear Dirac equation ( 12.21 ). 

Properties of solitary wave solutions 
Theorem 2.1. Let n > 1, f e C* 00 ^), /(0) = 0. 

1. Any solution (p u to J2.31 ) with some uj G (—m, m), is such that then for any /i < \Jm 2 — uj 2 one has W £ 
^(R",^), where (Q) = y/l + x 2 . 

2. For uj £ R\ [ — m, m], there are no solitary wave solutions satisfying (Q)^<j>u € L ca (W n ,C N ). 
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3. Assume that n < 3 and k G N. If n — 3, additionally assume that k = X. If f G C°° (R) satisfies 

f(v) = V k + o(ri k ), 

then there is Uq < m rac/i f/iere are solitary wave solutions (f> UJ (x)e^ lult to (12.21 > w/f/; w £ (woj 

This theorem is proved in Section[3] In Section [T2l we will also derive the nonrelativistic asymptotics of solitary 
waves as cj — > to; this will allow us to make conclusions on the spectra of linearizations at solitary waves. 

Properties of embedded eigenstates 

Consider the solution to (12.21 i in the form of the Ansatz ip(x, t) = (<j> u (x) + p(x, t))e~ lult , so that p(x, t) G C N is a 
small perturbation of the solitary wave. The linearization at the solitary wave 4> LLl (x)e~ lu ' t (the linearized equation on 
p) is given by 

d tP = JC(io)p, (2.4) 

where 

J=l/i, C(u) = D m - u) - MPMP - 2/ '{<f>tP<l>M« Re(C P ■ )• (2.5) 
Remark 2.1. The operator C(lu) is not C-linear because of the term with Re(0* p ■ ) (cf. Remark [L2l . 
Let O be an open subset of [— m, m] as in Assumption[T] 

Theorem 2.2. 1. If uj e O and X £ <J V (J C(uj)) fl iM, |A| < m + f/zen the corresponding eigenfunctions are 
exponentially decaying. 

2. There are no embedded eigenvalues beyond the embedded thresholds: 

a p {JC{u)) n i(R\[-m - \u\,m + = 0. 

Ifn = 1, then additionally ±i(m + ^ Up{J C(u))). 

We prove this theorem in Section 

Bifurcation of eigenvalues with nonzero real part 

By Weyl's theorem on the essential spectrum (see Lemma |4~TI below), the essential spectrum of JC(uj) is purely 
imaginary; the discrete spectrum of is much more delicate. Our aim in this paper is to investigate the presence of 
point eigenvalues with positive real part. The presence of such eigenvalues leads to the linear instability of a particular 
solitary wave: certain perturbations start growing exponentially. As uj changes, such eigenvalues can bifurcate from 
the point spectrum on the imaginary axis or - possibly - even from the essential spectrum of J L{iS). We will show 
that the bifurcations of point eigenvalues from the essential spectrum into the half-planes with Rc A ^ are only 
possible from the embedded eigenvalues, and that there are no bifurcations beyond the embedded threshold. There 
are only two exceptions, which could occur when the edges of the continuous spectra meet: bifurcations from A = 
at w = ±m (as in [CGG12]) or bifurcations from ±mi at uj = (as in [KS02]) do not have to start from embedded 
eigenvalues. 

Theorem 2.3. Let n > 1. Let (ujj)j^, ujj € O, be a Cauchy sequence, and let Xj 6 a p (/J C(ujj)). Denote 
ujf, := lirrij-i.oo ujj, and let Aj € CU {oo} be an accumulation point of '(Aj)jgN- Then: 

1. Xb 7^ oo. 

2. If Xb € iM, then \Xb\ < m + \oJb\- Ifn = L then moreover \Xb\ < m + \u>b\. 

3. If u>b € O and \Xb\ < m + \uJb\, then Xb G <J p (J 'C{uJb)), and there is a subsequence of eigenfunctions (Cj')j'eN 
corresponding to Xj G a p (^C(u>j)) which converges in L 2 to the eigenfunction £b corresponding to Xb G 
a p (J C(u>b))- 

Moreover, if there is a subsequence of (Aj)jgBj such that Xj Xb and Re Xj ^ 0, then 

(Cb,c(Lo b )Cb) = 0. 
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4. Define the zero order self-adjoint operator V{uj) : L 2 (R n ,C N ) -> L 2 (E", C N ) by 

V(uj) = C(uj) - D m + uj. 
IfuJb = ifi and there is s > 1/2 .smc/i that 

lim ||(Q) 2s ^(w J )|| i oo (K „ iEnd(C iv )) =0, 

r/ien A fc G {0; ±2mi}. 

If additionally Re Aj 7^ 0, fnen Aj, = ant/ moreover Xj = 0(m — \cjj\). 
This theorem will be proved in Sections|6]and|7] 

We are going to relate the families of eigenvalues of the linearized nonlinear Dirac equation bifurcating from A = 
with the eigenvalues of the linearized nonlinear Schrodinger equation. Let n < 3. Let u k (y), k G N (k = 1 if n = 3) 
be a strictly positive spherically symmetric function to the equation 

- ^u k = -^Au k -\u k \ 2k u k , u k eS{R n ), (2.6) 

so that u k (x)e~ lult with cj = — 77- is a solitary wave solution to the nonlinear Schrodinger equation 



2m 

1 



2m 

By (11.7b . the linearization at this solitary wave is given by 

9tP = jlp, 



where j 



1 

-1 



1+ 
l_ 



, and 



I- = A ~ A ~ uf, 1+ = A - A - (2* + l)u» (2.7) 
2m 2m 2m 2m 

where is a strictly positive spherically-symmetric solution to j2.6h 

Remark 2.2. By |BL83 , Example 1], equation j2.6i in K" with k G N has nontrivial solutions if and only if n < 3; if 
n = 3, one additionally requires fe = 1. 

Theorem 2.4. Lef n < 3, and let f G C°°(M), /(??) = r/ fe + 0(77*), itel If n = further assume that k = 1. Lef 
4> u e~ lut be solitary waves from Lemma [PI ( see below). Let ujj — > m, ana" /ef Aj € <J V {3 £.(ojj)), Xj — 0(m 2 — u> 2 ). 
Assume that the operator l_ introduced in (I2.71 l satisfies ) = {0}. 

7. AssMme f/iaf A = is an eigenvalue of )\ of algebraic multiplicity 2 + 2n (either n = 1, fc 7^ 2; or n = 2, 
A: > 1; or n = 3). 77ien Aj 7^ 0/or j e N implies that 

A h := lim -A — rr 7^ 0. 

j-s-oo m z — us* 

2. A b e (T p (jl) or A b G ct p (1_). 

3. If Tie Xj 7^ 0, j € N, fnen Af, new zero Krein signature: a corresponding eigenvector z G L 2 (M") satisfies 

(z,lz) = 0. 

Remark 2.3. By [CGNT08], the condition <7 p (l- ) = {0} is satisfied in the case n = 1, fe = 1 or k = 2, and also in 
the case n = 2, = 1. In the cases n — 1, fc > 3 or n = 2, > 2, or n = 3, fc = 1 the small solitary waves are 
already known to be unstable due to the presence of positive eigenvalue in the spectrum of the linearized equation; see 
HCGG12L 

Remark 2.4. The algebraic multiplicity of A = G cr p (jl) is at least 2n + 2. Indeed, since kerl_ = Span{ufe}, 
ker l + = Span{9jUfe ; 1 < j < n}, dimker jl = n+ 1. Moreover, each of these null eigenvectors, being orthogonal 
to the kernel of (jl)* = Ij, has an adjoint eigenvector. 
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Spectral stability of solitary waves of the nonlinear Dirac equation in ID 

We use the above results to prove the following spectral stability result. 

Theorem 2.5. Let f e C°°(R), /(O) = 0, /'(0) > 0. There is luq 6 (0,m) such that there is a family of solitary 
wave solutions <j) u] (x)e~ 1LJt to the nonlinear Dirac equation, 

id t ^ = D m ^ - f(i>*pi>)pi>, %l>(x, t) e C 2 , x e K, (2.8) 

and for each u> € (ujq, m) the corresponding solitary wave is spectrally stable. 

Above, D m = —iaJ^; + /3m, where a, j3 are self-adjoint matrices and satisfy a 2 = (3 2 = I2, {a, (3} = 0. 

Proof. We consider the family of solitary wave solutions w e~ lu * which is described in Lemma l3~5l below. 

Remark 2.5. In ID, this family is unique, in the sense that for each u> G (—to, to) \ {0} there is at most one solitary 
wave (see e.g. HBC12 I). modulo the translations in x and the U(l)-invariance. 

We assume that there is a family of eigenvalues Xj € cr p (JL(w)), with Rc Xj ^ 0. Then, by Theorem |2.2| Xj — » 
and moreover Aj := 3 — > At e cr(jl), where jl is the linearization of the cubic NLS in ID, ( dl.41 >) with n = 1 

and f{rj) = r\ + 0(77)). By Theorem |2.4| ([2l), A b := lim Wj ^. ro — gj— a belongs to the point spectrum of the cubic NLS 

linearized at a solitary wave; by ICGN T081 Fig. l,p = 3], this spectrum consists of A = only (note that there could 
be no embedded eigenvalues in ID due to the asymptotics of the Jost solutions), hence we must have A& = 0. 

By Theorem 12.41 rfTT). since the generalized null space of the linearization of the cubic NLS is four-dimensional, 
Aj = for all but finitely many j; thus, there is no sequence (Xj)j^ with the above properties. □ 

Remark 2.6. A slightly more careful analysis shows that small amplitude solitary wave solutions to the Dirac equation 
in ID with quintic nonlinearity and in 2D with cubic nonlinearity are spectrally stable. (In other cases, small amplitude 
solitary waves are linearly unstable HCGG12I .) 

Remark 2.7. While we try to exclude the bifurcations of nonzero-real-part eigenvalues from the essential spectrum, 
there is a possibility that, as uj changes, purely imaginary point eigenvalues bifurcate from the edges of the essential 
spectrum into the spectral gap (so that Re A = 0) even if when there are no embedded eigenvalues at the edges. This 
was noticed numerically in e.g. BBC 121 . Then a pair of purely imaginary eigenvalues could either collide at A = and 
turn into a pair of one positive and one negative eigenvalues (this collision is characterized by the Vakhitov-Kolokolov 
condition <9 W ||</> W || 2 = 0; see e.g. |Coml 1 1), or two pairs of purely imaginary eigenvalues could collide in the gap but 
away from A = 0, producing a two pairs of eigenvalues with nonzero real parts. 

3 Solitary wave solutions 
3.1 General properties 

Lemma 3.1. Let n > 1. Let ^ w e L 2 (M. n ,C N ), lu e ( —m, to) be a solution to ( 12.3I >. Then for any fi < \J m 2 — ui 2 
one has 

e^ Q) <j> u £ L°°(R n ,C N ). 

As before, (Q) is the operator of multiplication by yl + x 2 . 

Proof. For the sake of completeness, we choose to provide a proof of the above lemma. We will use the Combes- 
Thomas method, see [His00|. The solitary wave profile <fi u satisfies 

= D m ct> u - fttZ^Mu. (3.1) 

Pick /j, e (0, \/ m 2 — oj 2 ) and let e g (0, V m 2 — u> 2 — p). Due to the AssumptionQ] <f> u (x) are smooth functions of 
x, tending to as \x\ — > 00; so is f{4>* LU f34> u ). Therefore, we can write 

fifoPM = f (x,lu) + f 1 (x,lu), x e r, 
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where Fo and F\ are smooth and Fq has compact support, while sup^Rn weI l-fiO^kOI — e - F° r an y f G C 2 (R"), 
we have: 

e^(Z> m - w - Fifie-VeV^ = F Q pe^cj> u , 

or 

(-Dm - w - F x p + D oV >)e*ct> u = F /3e^ w . (3.2) 



We take <p(a;) = ^p(x), with p(x) = (x) := \/l + .t 2 , so that lim |Vp(cc)| = 1, ||Vp||z,~ = 1. By d3~2b . 

|sc|— >-oo 

e^i = (An - w - fiiS + uDoP^FoPeffa. 
The invertibility of £> m — a; — F\{5 + fiDop follows from the condition on /i and the estimate 

|| (An - w - Fij8)«|| > dist(w,±ro) - e, 

leading to 

||e^|| < (Vm 2 -u?-e- ^ph^'W F o 0e^\\ L oo UJ = l^feML . 

V m — or — e — fi 

Since e Aip ^ u satisfies the elliptic equation d3.21 i, boundedness of e MP u in L 2 (R™, C^) leads to its boundedness in 
H k (R n , C N ), for any k e N. By the Sobolev embedding theorem, ||e w ^ w ||i,»o is bounded. □ 

For solitary waves associated to embedded eigenvalues the situation is, in many respects, different. We can first 
notice that with Hardy type estimates, we have 

Lemma 3.2. Let n > 1. Let w£l \ [—to, to]. If a solitary wave <j) bJ (x)e~ lu}t associated to uj is such that 

(Q) 1/2+£ <^ G L°°(R n ,C N ) 
for some e > 0, then W decays faster than algebraically in x. That is, for any K > 0, 

(Q) K ^eL°°(R n ,C N ). 

Proof. Since the nonlinearity / <E C°°(K) from dZ2l ) satisfies /(0) = 0, there is g e C°°(R) such that /(ry) = 775(77). 
Then, from Lemma lA31 for any s > —1/2 there is a constant c such that 



UQVcbjm < c\\(Qr +i (D m - a)<m < c\\g{rpmQ) 1/2+e rP(Q) 1/2+ ^mQ) s - 2e ^i 

which bootstraps immediately. □ 
The situation is even more dramatic if one applies the Carleman-Berthier-Georgescu inequalities from AppendixlBl 

Lemma 3.3. There are no solitary wave solutions (f> u (x)e~ with u> G M \ [— to, to] smc/z that 

(Q) 1/4 <^ GL°°(M",C JV ). (3.3) 

Proof. The proof goes in three steps. 

1. Let w £ R, |w| > to. Denote D a = {x € l n ; |x| > i?}. We have the following: There exists R a > 0, t > 
and C(Rq) > such that for any i? > i?o, for any it € Hq{Qr, C n ), and for any r > to there is the inequality 

l|e T|Ql u|| i2 (^c«) < ^^Hr^e-IOI^-/^*^) - w)u||i»(fl«,c")- 
Indeed, by Lemma |B~4l there are C, R < 00 such that for any u € H^(f2ji, C N ) and r > 1 

r|| e " u || < C||r 1 / 2 e Tr (I? m -A)u||, 
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and thus 

rlle^lU^c*) < 9.\\ r ^ e rr {Dm _ f^J^) _ u)u\\ L . {QRtCN) + ^||e^V(C/^c>IU^ H ,c*)- 
So that from the assumption, we obtain the claim for sufficiently large r. 

2. Now we can prove that a solitary wave </> w with \uj\ > m which satisfies ( 13.31 ) is smooth and with compact support. 
Indeed, the smoothness follows from AssumptionQ] Then consider Vj = lytf) with r/j := rj(-/j), where 77 is smooth 
and satisfies < rj < 1, identically equals 1 outside the ball of radius 2, and identically equals zero inside the ball of 
radius 1. From the previous step, it follows that, for r > 1 and for j > i?o sufficiently large, 

l|e T|Q mi i2( ^,c«) < ^^|||Q| 1/2 e T|QI (A„ - f(KPM-^)vj\\^n R ,c^y 
Therefore, ||e T l«lt; J -||i a(fljl ,c*) < 22 ^ sl || \Q\ 1/2 e T ^a • (Vt^Jl^x"), which implies that 
e T3k Hu\\mn 3k ,c») < const — )^- L e r2k \\<l> u \\ t nn k nB ak ,c"), 

where 

={i£ M™ ; |x| < R}. 

Since t could be arbitrarily large, one concludes that W is identically zero outside of the ball B^k- 

3. Now we need the unique continuation principle for the Dirac operator HBG87I . 

Lemma 3.4. Let 57 an open connected subset of "R" . Let v : fl — > M + be in Lf oc (fl). Ifip G _ff^ c (f2) satisfies 

\(D ip)(x)\ < v(x)\ip(x)\ a.e. in fl 

and ip{x) = in an open non-empty subset ofVt, then ip = 0. 

We refer to [BG87 Appendix] for the proof. Although it is written in dimension n — 3, the key part [BG87 
Appendix, Lemma 1] is true in any dimension. 

From Lemma (13.41 1. the solitary wave ifi^ is identically zero. □ 

Remark 3.1. The statements of the above lemmas can be improved with the stronger assumptions on /. For instance, 
if f(v) — 0(\rj\ k ), then there are no solitary wave solutions (f> U j(x)e~ lu ' t with u> 6 K \ [—m, m] such that 

(Q) 1/2k (t>u e L°°(R n ,C N ). 

Lemmas |3 . 1 1 and l3~3l complete the proof of Theorem l2.ll 



3.2 Solitary waves in the nonrelativistic limit 

In this section, we will study asymptotics of solitary waves in the nonrelativistic limit lu — s- m. The proof of the 
following lemma is given in [CGG12]. 

Lemma 3.5. Let n < 3. Let 

/eC°°(R), f(v) = v k + o(v k ), fceN. (3.4) 

Ifn — 3, additionally assume that k = 1. There is luq < m, dependent on n and f, such that for lu G (luo, tn) there 
are solitary wave solutions 4> bJ {x)e~ luJt , with W satisfying 

lu^ = D m (j) - f((f>Zf3<f>u)P4>u- (3.5) 
Moreover, introducing the projections onto the "particle" and "antiparticle" components, 

n P = ^{l + (3), 77 A = i(l-/3); cf> P {x) = n P <t>(x), 4>a{x) = n A 4>(x), 
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we have for e = %/ m? — oj 2 

4> P {x) =ei& P (ex)+0 H 2(e 2+ i), </> A (x) = e 1+ ^S A (ex) + O m (e 3+ ^ 
where <Pp{y), ^aIv) satisfy 

1 S P = -J—A&p - \Sp\ 2 kS P} & A — -^—(—ia ■ V y )Sp. 



2m 



2m 



2m 



One can choose 



(3.6) 

(3.7) 
(3.8) 



Remark 3.2. As we mentioned in Remark 1221 equation ( 12.61 ) in K™ with k € N has nontrivial solutions if and only if 
n < 3; if n — 3, one additionally requires k = 1. Absence of corresponding solutions to ( 13.7b does not allow us to 
construct small amplitude solitary wave solutions (in the nonrelativistic limit oj m) to the nonlinear Dirac equation 
in K 3 with k > 2 and in R", n > 3, with any k e N. 



# P (y) = nu k (y), 

where n G C^, ||n|| = 1, ant/ Ufc S <S(K") is a strictly positive spherically symmetric solution to 



4 Linearization at a solitary wave 

4.1 Essential spectrum of the linearization operator 

Denote 

C-(u)=D m -u-f(<l>lP<l> u )p, % :=V(£-(lo)) = H 1 (W\C n ). 

Introduce 



a J = 



i 1m a J i Re a? 
— iRea^ ilma^ 



Let 



Denote 



1 < J < n; 



P = 



Re <^ (x) 
Im (x) 



Re ^ - Im 
Jmj3 Re/3 



J 



In 
-In 



(4.1) 



(4.2) 



L_(w)=D m -a;-/(c|>iP<|>Jp ) L(w) = L_(w) - 2/'(cp* pcp w )(cb* (3 • )|3cb w . (4.3) 
Above, D m = — ia • V + mp. If <f> u (x)e~ lut is a solitary wave solution to (12.2b . then (p w satisfies JI (w)cp w = 0. 



The linearization at the solitary wave ( 12.4b takes the form 

<9 f p = JL(w)p, p(ar,t) = 



Re j o(x,i) 
Im p(x, t) 



»2N 



(4.4) 



Both I and L act on H 



lflB>n to 2 AT \ 



by C-linearity, we extend them onto X = H ( 



p2iV\ 



Theorem 4.1 (Weyl's essential spectrum theorem, |RS78|, Theorem XIII. 14, Corollary 2). Let A be a self-adjoint 
operator and let C be a relatively compact perturbation of A. Then: 



• B = A + C defined with D(B) = D(A) is a closed operator. 

• <T CSS {B) = <T oss (A). 



Applying Theorem l4. 1 I to A = iJ(D TO — oj), B = i(JL(w) — A), we conclude: 
Lemma 4.1. a css (]L(oj)) = i(M\(-m + \u\,m - |w|)). 
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4.2 Independence on the choice of Dirac matrices 

Let P € End(C Ar ) be a hermitian projector which commutes with all Dirac matrices: P 2 = P, P* = P, [P, a 3 ] 

TRe P — 1m P 

0, 1 < j < n, [P,f3] = 0. Denote JT = P3£, %i = (1 - P)3C . Denote P 



ImP ReP 

with <x> and with |3. Since P is C-linear, P commutes with J. The relation P<j> = <f> leads to P(p = cp. 



P commutes 



Lemma 4.2. If ' P<p u = 4> u , then the linearization at cf> u e lu>t satisfies: 

1. (7 ess (JL) = a 0SS (JL_) = a css (JI-lx ) = cr es8 (JL| Xi ) = i(M \ (-m + m - 

2. a(JL| Xi ) = a(JL_)ciK; 

3. a p QL\ x ) = CTp (JL_) U a p (JL| Xo ), a p (JL| x ) \ iR = a p (]L\ Xo ). 

Proof. The statement about the essential spectrum is immediate due to Lemma |4~T1 The inclusion cr(JI ) C iR 

follows from I being self-adjoint (with its spectrum a subset of R) and commuting with J (which has ±i as its 

eigenvalues when acting on C 2N = C <E>m R 2N ). 

JL acts invariantly in both Xo and Xi. Due to this invariance, one has cr(JL) = <r(JL| x ) U er(JL| x ). Moreover, 
JL| Xi = JL_| Xi ,hencecr(JL| Xi ) = cr(JL_| Xi ) C iR. We conclude that cr(JL) \ iR = cr(JL| Xo ) \ iR. ' □ 

Remark 4.1. Since cr( JL_) C iR, it follows that er( JL| Rangc p ) \ iR = a{ JL) \ iR. Therefore, the linearization of 
( 12. 2t at the solitary wave <p UJ (x)e^ lult has the same point spectrum away from the imaginary axis as the linearization 
at its embedding into the space of spinors of higher dimension. In particular, if there were a family of eigenvalues 
Xj of JL(wj)| Rallgo p bifurcating from X b G (T css (JL(a;&)| RanBC p ) such that X b <£ cr p (JL| Rangc p ), then there would be 
the same family of eigenvalues of JL(u;) bifurcating from € (T css (JL(ojb)), X b ^ cr p (JL(u;{,)). Therefore, before 
studying the bifurcations of point eigenvalues, we can first embed the Dirac equation and a particular solitary wave 
solution into the spinor space of higher dimension. 

Let us show that there is no dependence on which embedding we choose, as long as n is odd. 

Lemma 4.3 (Dirac-Pauli theorem). Let {a J , 1 < j < n; /?} and {a , 1 < j < n; /3}, be two sets of the Dirac 
matrices of the same dimension N: 

{a\a k } = 26 jk , {a J ,/3} = 0; {&>, a fc } = 25 jk , {&J}=0. 

1. Let n = 2d + 1, d G N. There is an invertible matrix S such that 

& = S"VS, 1 < j < n; {3 = S _1 /3S. 

2. Let n = 2d, d G N. There is an invertible matrix S and a G {±1} such that 

& = aS- l oPS, l<j<n; = aS~ l 0S. 

Above, one could choose S to be unitary. See HPau36l lvdW32| [Dir2 8 1 , |Tha92 Lemma 2.25], and also HKes61| 
Theorem 7] for general version in odd spatial dimensions. 
Let us give the sketch of the construction from [Fed96|. 

Proof. Let us remind the standard construction of the irreducible representation of the Clifford algebra with 2d gener- 
ators. Let ej and fj, 1 < j < d, be the generators of the Clifford algebra Chd- 

{e j ,e k } = {f j ,f k } = 26 jk , {e 3 ,f k } = 0; 1 < j, k < d. 

Define Zj = |(ej + ifj), z* = |(-ej + ifj). Then z 2 = (z*) 2 = 0, {zj,z* k } = S jk . The operators Zj and z* are 

often referred to as the operators of creation and annihilation. Define the "vacuum vector" po = Ylj=i ( z j z j ) e Chd, 
and let S = {vpo ; v 6 C^d} C C^d be the left ideal of pq. Since Zjpo — for 1 < j < d, the elements of the 
form (z\ ) ai . . . (z*[) ad po, with oi, . . . , dd € {0, 1}, form the basis of the spinor space S; hence, dim S = 2 d . The 
action of Cl-2d m S is the only irreducible representation of C^d- Any other representation of Cl^d is isomorphic to 
this representation or to its several copies. 
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In the case of the even-dimensional space, x e M™ with n = 2d, we need consider the Clifford algebra C^d+i- 
We will denote its generators by ey, 1 < j < d + 1, and fk, 1 < k < d: 

{e j ,e j ,} = 25 jf , {fk,fk>} = 2 <W, {e,-,/ fe } = 0; 1 < j, j' < d+ 1, 1 < k, k' < d. 

We consider C^d+i as embedded into Cl2d+2, adding one more generator fd+i which anticommutes with ej, 1 < 
j < d+l and with fj, 1 < j < n, and also satisfies = 1. We introduce = ~(ej+i/j) andzj = |(— ej+ifj), 

1 <i < d+l, and build the representation space S whichis the left ideal of the element z j( z j)*> dim S = 2 d+1 . 
The action of Cl2d+i in S is reducible. We decompose S into S = So ® Si, where So consists of elements which 
can be represented as the products of a certain number of ey, 1 < j < d + 1, and / 3 , 1 < j < d, while Si consists 
of elements which can be represented as the products of some e 3 , 1 < j < d+l, some fj, 1 < j < d, and one 
copy of the element fd+i- Then Chd+i acts invariantly in So and Si. Moreover, each of these representations is 
irreducible since dim So = dim Si = 2 d , while the only irreducible representation of the subalgebra Clid C Clid+i 
is of dimension 2 d . The representations of C^d+i in So and Si are not isomorphic. □ 

Remark 4.2. The representations of C?2d+i in So and Si are related as follows. The left multiplication by fd+i 
defines the mapping 

is : So ->• Si, pi : g ^ fd+iQ € Si, 

with the inverse 

pT 1 : Si -> S , M 1 : h ^ fd+ih G So- 
Due to the relations 

fd+iejfd+i = -ej, 1 < j < d + l, fd+if 3 fd+i = -fj, 1 < i < d, 
the action of e^, 1 < j < d + 1, and fj, 1 < j < d, by the left multiplication on So and on Si are related by 

V° Ls (ej)o [j,- 1 = -L Sl (ej), l<j<d+l, fx o L So (fj) o (J.' 1 = -L Sl (fj), 1 < j < d. 
Lemma 4.4. The spectrum of JL is symmetric with respect to the real and imaginary axes. 



Proof. Due to Lemma I47T1 it is enough to consider the discrete spectrum. The inclusion — A € er(JL) follows from 
(— JL)* being the conjugate to JL: 

(-JL)*=L*(-J)*=LJ = J- 1 (JL)J. 

Let r be a self-adjoint matrix which satisfies T 2 — 1, anticommutes with a J , f3, and also satisfies T<f> — <fi. (Such 
a matrix exists if we embed the Dirac equation into the spinor space of sufficiently large dimension; by Lemma l4~2l 

["ReT -imr" 



this does not change the spectrum.) Then T = 



is a self-adjoint matrix which satisfies T 2 = 1, 



anticommutes with a 3 , (3, commutes with J (since T is C-invariant), and also satisfies Tcp = cp, where cp = 



Re< 
Pm ( 



Then rjIX 1 = -JL, showing that A e cr p (JL) implies -A E cr p (JL). □ 

4.3 Point spectrum of the linearization operator 

Lemma 4.5. Let a be an hermitian matrix anticommuting with a 3 , 1 < j < n, and with j3. Then a°(j) is an 
eigenfunction of C- and of C, corresponding to the eigenvalue A = — 2d. 



Remark 4.3. If n = 3, one can take a — a 1 a 2 a 3 /?; with the standard choice of Dirac matrices, a = i 



-h 
h 



Proof. Since a° anticommutes with a 3 , 1 < j < n, and with (3, and taking into account ( 13.11 ), we have: 

C-a°4> = (D m -w- f((j)*/3<j))f3)a (j) = a°(-D m - lu + f{<j>* (3<j))0)<j) = - 2uj)(j) = -2ua°<j). 



Since a and /3 are Hermitian, 2Re[0*/3a°0] = + </>*/3aV = 0*{/3,a°}0 = 0; therefore, one also has 

= -2ua°d>. □ 
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It follows that the linearization operator has an eigenvalue 2wi: 

2ui G a p (JC(uj)). 

Since cr(JL(w)) is symmetric with respect to R and iM, for any f(rf) in d2.2| ) and in any dimension n > 1, we arrive at 
the following result: 

Lemma 4.6. ±2u;i e er p (JL(u>)). 

Remark 4.4. For > m/3, the eigenvalues ±2wi are embedded into the essential spectrum. 

Remark 4.5. The result of Lemma l4~6l takes place for any nonlinearity f (ip* (3ijj) and in any dimension. The spatial 
dimension n and the number of components of ip could be such that there is no matrix a which anticommutes with 
Oi 3 i 1 < J < n>, and with j3; then the eigenvector corresponding to ±2wi can be constructed either using the spatial 
reflections. Alternatively one can double the size of the spinors so that there is a desired matrix a ; by the results of 
this section, this does not change the spectrum of the linearized operator. 

5 Properties of embedded eigenstates 

5.1 Decay of embedded eigenstates before the embedded threshold 

Lemma 5.1. Let A € ct p (JL(oj)) satisfy A £ iM, m — \u>\ < |A| < m + |cj|. Then the corresponding eigenfunction is 
exponentially decaying. 

Remark 5.1. The proof for the case n = 1 follows from the properties of the Jost solutions (any eigenfunction 
can be decomposed at x — ±oo over the exponentially decaying Jost solutions, and only the ones exponentially 
decreasing at ±oo can participate; thus, the exponential decay of Jost solutions translates into the exponential decay 
of eigenfunctions). 

Proof. The eigenfunction I e L 2 (W\C 2N ), \\C\\ = 1, which corresponds to A, satisfies (D m - u + AJ)£ = -VC 
Applying II*, we get: 

(D m -u + iA)C+ = -I7+VC, (D m -u- iA)C = -ITVC 

Assume that iA is of the same sign as ui (the other case is treated verbatim by exchanging the treatment of t^); then 
uj — iA is in the spectral gap of D m , hence we can write 

C+ = -(D m -w + iA)- 1 n+VC 

The Combes-Thomas method |HisOO|, as in the proof of Lemma |3.1| shows that £ + is exponentially decaying. 

Now we turn to C . Due to the exponential decay of W ( x) (see Lemma [3TTb and hence of V( x, uj), we can choose 
t > small enough so that 

sup ||r 1/2 e T|:r| n"V(.T,a;)|| End(C 2iY ) < oo. 

Let 0(x) E C°°(IR n ) be such that = 1 for \x\ > 2, = for |ar| < 1. Applying Lemma|H31 we have: 

||e Tr e(x)^|| < C||r 1 / 2 e rr (D m -a;-iA)6)(.x)^|| < C\\0{xy/ 2 e Tr U- VC|| + C\\r^ 2 e Tr [0(x), D ]C|| • 

It follows that e^l^l C is bounded in L 2 . 

Taking derivatives in x and using theboundednessof Cin7J s (R",C 2Ar ),Vs G N, we conclude that e T W0(x)Z~ (x) £ 
H s (WL n ,C 2N ), Vs G N. The exponential decay of £T and hence of C follows. □ 
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5.2 Absence of embedded eigenvalues beyond the embedded threshold 

Lemma 5.2. Assume that there are e > and C < oo such that |V(x)| < Ce _e l x l, x G R™. If A G Op(JL) n iR ami 
| A| > m + then for any N > fnere is Cn > swc/z f/iaf fne corresponding eigenfunction C satisfies 

\Z(x)\ < C N e- Nlxl , xeR n . 

Remark 5.2. The proof for the case n = 1 follows from the properties of the Jost solutions: while any eigenfunction 
can be decomposed at x = ±oo over the exponentially decaying Jost solutions, there are no decaying Jost solutions at 
thresholds A = ±i(m + |w|). 

Proof. The proof is a bootstrap argument based on Lemma lB~4l As both uj ± Im A have modulus bigger than m, from 
this lemma we deduce that if e Tr (j(D m — uj) — \)r)C, is square integrable for any bounded smooth function r\ with 
support outside some sufficiently large ball then e Tr r)t, is as well. As 

(J(D m - w) - A) ( V Q = Jia • (Vtj)C - JV^rtf, 

the bootstrap starts with the exponential decay of V. □ 

Lemma 5.3. Let n > 1. Fix uj G O. There are C < oo and Rq > 0, dependent on n, uj, and A, such that for any 
R> Ro and any u G Hq(S1r, C n ), 

T 1/2 \\e Tr u\\ < C\\r 1/2 e Tr (C(uj) ± iA)«||, r > 1, 

where £(uj)u = D m u -urn- fifcPMu - 2/'(>* (3^) Re(0* f3u)<t> u . 

Proof. This is the adaptation of the Carleman-Berthier-Georgescu estimates. By Lemma [B~4l for r > 2m and any 

u € H^(f2 R ,C N ),one has: 

2 

|| ((w T iA) 2 -m 2 + T —) 1/2 e rr u\\ < \\p(r)e Tr (D m -u± i\)u\\ 

< || A /(r)e^( J C(^) ± iA)w|| + ||/4r) e -(/(0*/3^) U + 2/'(0*/30) Re(^/3 U )^>)|| , 

with 

, . (n (ui^Wfr „ \ 1/2 

Note that (uj =F iA) 2 — m 2 > since A € iR, | A| > m+ \uj\. Therefore, once u> and A are fixed, if R > is sufficiently 
large, there is C < 00 so that 

\\e rr u\\ <CMr)e Tr (£(uj)Ti\)u\\, uGH^(Q R ,C N ), t > 1. 

□ 

Lemma 5.4. Let n > 1. 77ie operator JL(w) /ias no embedded eigenvalues A 6 iR w/fn |A| > to + 

Remark 5.3. For the case n = 1, the proof follows from the analysis of the Jost solutions, and the conclusion is 
stronger: the operator JL(w) has no embedded eigenvalues A G iR with |A| > to + \u>\. 

Proof. Assume that A G iR, |A| > m + \uj\, is an embedded eigenvalue of JL(w), with C G L 2 (R n ,C 2N ) the 
corresponding eigenvector: 

- JAC = D m C - w£ - /((p*(3(p)(3C - 2/'(cb*(3cp) Rc((p*|3C)|3cp. (5.1) 

Let 6* G C°° (R™) be a smooth radially symmetric cut-off function with support in the closure of f2 R .+i an d with value 
1 in By Lemma 1331 

Vr>l, \\e Tr 0C\\ < -^||r 1/2 e rr (JL- A)6>C||. (5.2) 
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Since JU = AC we have (JL - X)OC = [JL, 6>]C By IrO) . 

Vr>l, Vfc>i?, ||e Tr 0C|| < -^||r 1/2 e Tr [JL,e]C||. 
Taking into account that d Xj and hence [JD, 0} are zero outside of the ball we conclude that 

Vr> 1, Vfc >fl, He^ClU^^^) < ^ c (-B+a)T|| r i/2 V ©|U-I|C|U a (B* +a ,c^)- 

Since r > 1 could be arbitrarily large, we conclude that supp £ n Dr+2 — 0- The unique continuation principle, 
Lemma [3741 ensures that £ = 0, contradicting our assumption that there were an embedded eigenvalue A £ iR, 
|A| > m + |w|. □ 

Theorem |2.2| follows from Lemmas |5 . 1 1 and [5~4l 



6 Bifurcations of eigenvalues from the essential spectrum 

Lemma 6.1. Let J £ End(C 2Ar ) be skew-adjoint and invertible and let L be self-adjoint on L 2 (R",C 2JV ), If X £ 
a p ( JL) \ iR w/f/i the corresponding eigenvector C f/ze« 

(W) = o, (CK) = o. 

Proof One has JL£ = AC L£ = AJ -1 C hence 

(cu) = A(cr 1 c. (6.i) 

Since (C LC) € R and (C J X C) G iR, the condition Re A ^ implies that both sides in ( 16.11 ) are equal to zero. □ 

Remark 6.1. If an eigenvector C corresponding to A £ cr p (JL) satisfies (C LC) = 0, we will say that A has zero Krein 
signature. The Krein signature is only interesting for A € iR since, according to Lemma |6*T1 all eigenvalues of JL 
with nonzero real part have zero Krein signature. 

We will formulate our results for the operator L(cj) = D m — u + V(u/), uj £ [— m, to], with V operator-valued, 
with V(w) zero order and self-adjoint, for each uj £ [—to, to]. Note that L(cj) is not necessarily a linearization at a 
solitary wave of the nonlinear Dirac equation. 

Lemma 6.2. Let n > 1. Let J £ End(C 2Ar ) be skew-adjoint and invertible and let 

L(w) = D„, — uj + V(cj), uj £ [—to, to], 

with V : [— m,m] — > L°°(R™, End(C 2Ar )) a zero-order operator-valued function which is self-adjoint for each 
uj £ [—to, to]. Lef ujb £ [—to, to], one/ assume that there is e > smc/i f/iaf 



limsup ||(Q) 1+e V(w)|| L <» (Rn!End( c2iv )) < oo. (6.2) 



a+£ 

,End(C 2JV )) 
a;— >a;b 



Letujj £ (—to, to), cjj — )> Wfc. Let Xj £ er^JL^)) be a family of eigenvalues such that T\e Xj ^ 0, Aj — ^ A/, £ iR. 

j-s-oo " ' ' J-i-oo 

G {0; ±to}, additionally assume that 

A fc ^±i(TO-|w fc |). (6.3) 

Then 

X b £ a p (]L(uj b )). 

Moreover, there is an infinite subsequence of eigenfunctions (Cj)jen corresponding to Xj £ er p (JL(wj)) which con- 
verges to the eigenfunction corresponding to Xb, and this eigenfunction Ci satisfies 

(C,lmc) = o, (c,JO=o. 
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Remark 6.2. In this lemma, we do not include bifurcations from thresholds in full generality, but if we replace ( 16.21 ) 
by a stronger condition that there is e > such that 

limsup |KQ) 2+E V(w)|| i oc (R „ iEnd(C 2jv )) < oo, (6.4) 
then the lemma remains true without any restriction on cut, or Ah- 

Remark 6.3. The conclusion of the lemma is trivial for Ab G iffi. with |A&| < m — \tOb\, when Ab is not in the essential 
spectrum, and the inclusion A/, G <7 p (JL(u;&)) follows from the continuous dependence of eigenvalues in the discrete 
spectrum on a parameter, to. 

Remark 6.4. Due to the exponential decay of U proved in Lemma [3~T1 the assumption ( 16.21 ) is trivially satisfied for 
linearizations at solitary waves with Wf, ^ dO. Due to ( 13.61 . the assumption ( 16.21 ) is also satisfied for the linearizations 
at solitary waves of the nonlinear Dirac equation from Lemma 13.51 for ojf, = m. On the other hand, according to 
asymptotics ( 13.61 ). one sees that the condition ( 16.4b is not satisfied for the linearizations as to — > m. The case u>b = m 
is studied in detail in Sections|7]and[8] 

Remark 6.5. By Lemma 1631 if Ah = ±i(m + \tOb\) and w& G" {0; ±m}, then one must have Ah G <7 p (JL(wb)). This 
excludes bifurcations from ±i(m + \to\), \u)\ < m, in the case n — 1, since, as the analysis of the Jost solutions shows, 
there can be no embedded eigenvalues at ±i(m + |w|). 

Proof. Let (Cj )j*eN t> e a sequence of unit eigenvectors associated with eigenvalues Xj, so that ]J.{uij)C,j = XjCj. It 
follows that 

(D m - Wj - = -V(iJj)Crj. 

Let II ± = |(1 =p ij) be the projectors onto eigenspaces of J corresponding to ±i. We denote Cf = Tl ± Cj- By 
Lemma|63l0= (CjJCj) = i\\tf\\ 2 - i\\ZJ\\ 2 , j € N, while 1 = ||^|| 2 = ||£t|| 2 + ||CJ|| 2 , j g N; we conclude that 
11^=11 =1/2 and 

((D m - Uj ) T = -n ± V(o; J -)Cr 

If the condition ( 16.3b is satisfied, then either iAh — Wb or — iAh — uib is not a threshold of D m . Let us assume that 
iAh + u)b is not: 

iA b + w h 7^ ±m. 

Then iAj + uij ^ ±m if j is large enough. From Lemma !a31 we deduce that for any s > —1/2 

II (Q)% 3 \\m < 4 (Q) s+1 ((D m - Uj ) - iA^tfW = C\\ (Q> s+1 n+V(c J )C, II , 
for some C < oo which does not depend on j. Hence, for s — e, we have 

\\(Q) s Cf\\m<cC, j'eN. 

It follows that is precompact in L 2 (R™, C 2JV ), and we can choose a subsequence which converges to a vector 

of norm lim^oo = 1/2. As (CJ) is weakly convergent, we conclude that there exists a nonzero weak limit 

which is necessarily an eigenvector. □ 

Lemma 6.3. Letn > 1. Let to j G (— m,m), LUj — > ujb G [— rn, m]. Let Xj G (T ( j(JL(w :; )) be a family of eigenvalues 

j->oo 

swc/z f/iaf Aj — >■ Ah G iM. Assume f/iaf Ab is not at any of the thresholds, 

I A & I 7^ m±ujb, 

and that there is e > smc/; that 

limsup||(Q) 1+£ V(a;)|| L oo (R „ iEnd(C 2jv )) < oo. (6.5) 

Then Ah G cr p (JL(a;h)) and there is a sequence of eigenfunctions Cj corresponding to Xj G <7 p (JL(a/j)) which 
converges in L 2 to the eigenfunction £ b corresponding to Ah G a p {J\-(uib)). 

If moreover, Re Aj 7^ 0, then the corresponding eigenfunction Cb satisfies (£ b , L(wh)Ch) — 0, (tf,, JCfc) = 0- 

Proof. If iA — cj and — iA — cj are not thresholds of D m then the proof of Lemma |6T2l shows the relative compactness 
of both (0^~)jeN and (4>j)j e ?%. The rest is proved exactly as in Lemma |6\2l □ 
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7 Bifurcations from the essential spectrum of the free Dirac equation 

The limiting absorption principle for the Dirac operator has been established in [ Yam73 ] . We need a slightly stronger 
version for A G JC where JC is a closed set which does not contain the thresholds ±to, but is not necessarily compact: 

Lemma 7.1. Let JC c C be a closed set such that ±m G" JC. Then for any s > 1/2 there is C < oo such that 

II«IUl.(b»,c*) <C||(An-A)u|U i(K n >CW) , VAg/C, Vu € Hl(R n ,C N ). (7.1) 
Proof. For /C compact and A G /C \ R, this is the result of [ Yam73). For |A| > m + 1, 1m A 0, this follows from 

(D m - A)- X u = {D m + A)(-A + m 2 - A 2 ) _1 m, 

where one has 

||(-A + to 2 - A 2 n| L? ^ s = 0(|A 2 - m 2 |" 1/2 ) 
by | |Agm75| Remark 2 in Appendix A], and 

\\Dm + M\H*.^Hi. < C(l + |A|). 

This proves that for some C < oo which depends on to > 0, s > 1/2, and JC C C but not on A, one has 

\\{D m -X)- 1 v\\ Hl _ s <C\\v\\ L 2, \eJC\R, veL 2 s , 

yielding CO} for Jm A ^ and for any u G H] (R n , C N ). By continuity, O also holds for A e JC. □ 

Let us consider families of eigenvalues in the limit of small amplitude solitary waves, which may be present in the 
spectrum up to the border of existence of solitary waves: u> — » ujb E {±m}. This situation could be considered as the 
bifurcation of eigenvalues from the continuous spectrum of the free Dirac equation. 

Let J g End(C 2JV ) be skew-adjoint and invertible, and let 

L(cj) = D m — lu + V(x, oS), u) £ [— m, m], 

with V : [—to, m } -> L°°(R",End(C 2Ar )) an operator-valued function of cj e [— m, m], with V(w) a zero-order 
operator. 

Lemma 7.2. Lef uif, = ±m. Assume tliat there is s > 1/2 such tlzat 

limsup||<Q) 2 'V(w)|| L =. ( R»,End(C"')) =0- (7.2) 

Let ojj G O, uij — > Wh. Tjf Aj G Op(JL(cjj)), tlien the only possible accumulation points of {Xj ; j G N} are 
A = {0; ±2mi}. 

Remark 7.1. In this lemma, V(w) is not necessarily self-adjoint. 

Remark 7.2. By (13. 6K the condition d7.2t is satisfied for solitary waves from Lemma |331 in the nonrelativistic limit 
uj — > m, with any s < 1 

Proof. Let /C C C be a closed set such that ±m </ /C. According to Lemma ITTl there is the limiting absorption 
principle for the free Dirac operator D m = —ia ■ V + (3m, so that the following action of its resolvent is uniformly 
bounded for z G /C\M: 

(Dm-z)- 1 : L 2 s (R n ,C N ) — > L 2 _ s (R n ,C N ), s > 1/2, z G /C\M. (7.3) 

Now let V C C be an arbitrary closed set which does not contain ±2(to ± Wf,)i. To prove the theorem, we need to 
show that for to sufficiently close to lo^ there is no point spectrum of JL(o>) in V. Let to be close enough to uj^ so that 
V does not contain ±i(m ± uj). One has limui^oo L(cj) = D,„ — to. Since the eigenvalues of J are ±i, the operator 
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J(D m — to) can be represented as the direct sum of operators \(D m — to) and — i(D m — u). By ( 17.3b . the following 
map is bounded uniformly for z G V\iK: 



(](U m -uj)- z) : L 2 (R n ,C 2N ) ^ L 2 _ s (R n ,C 2N ), s > 1/2, z e V\iK. (7.4) 
The resolvent of JL(w) is expressed as 



(JL(w) - z) -1 = (J(D ro -w)-z) 1 —j. (7.5) 

1 + JV(j(D m - u) - z) 



Thus, the action (JL(w) - z)" 1 : L 2 (R n , C 2iv ) -> L?. a (R n , C 2N ) is uniformly bounded in z G V\iE as long as the 
operator V(w) : L^ S (M™, C 2Ar ) -> L 2 (]R", C 2Ar ) of multiplication by V(x, w) has a sufficiently small norm; it is 
enough to have 

l|V|| i ^ L j||(j(D rn -a;)-z)" 1 || L ^ i 2_ s < 1/2. (7.6) 
Due to the bound on the action J7.41 , the inequality ( 17.6b is satisfied since 

J^S ll^( w )ll i - s ( R "^ 2JV )->L2( R „ iC 2N ) < Jim ||(Q) 2;i V(a;)|| i oo( R , l!End ( C 2jv)- 1 = 

by the assumption of the theorem. □ 
Lemma 7.3. Assume that there is s > 1/2 such that 

J™ ll(Q) 2s V(w)|| L =e (Rn! End(c 2 «)) = 0. 

Tfwj — > lui, € {±m}, Aj G <7p(JL(cdj)), ReAj ^ 0, f/;en lim^oo A 3 = 0. 

Proof. Since JL(cj;,) = J(D m — cji,), is a differential operator with constant coefficients and with nondegenerate 
principal symbol, it has no point spectrum. Then Lemma loT2l provides the conclusion. □ 

Remark 7.3. Thus, bifurcations at cj/, = ±m off the imaginary axis are not allowed. At the same time, bifurcations at 
lji, = ±m from Xb = ±2mi are possible and indeed take place for the linearizations at solitary waves: by Lemma l4~6l 
there are families of eigenvalues A = ±2wi G cr p (JL(w)). 

Let J G End(C 2Ar ) be skew-adjoint and invertible, and let 

L(w) = D m — us + V(x, lo), uj G [— m, m], 

with V : [—to, to] — > L°°(R n , End(C 2JV )) an operator-valued function of u> G [—to, to], with V(cj) a zero-order 
operator which is self-adjoint for each cj G [—to, to]. 

Lemma 7.4. Let ujb = ±to ant/ assume that there is C < oo smc/i f/iaf 

||V(w)||i«,( R n )Elld (c2W)) < C\uj- uib], wG[-m,m]. (7.7) 

Lef G (—to, to), j G N; cjj — > cj&. If there is a sequence Xj G u(JL(cjj)), 5mc/; f/;af Re Aj ^ 0, lim Aj = 0, then 

j-tao 

Xj = 0(\u)j -u b \). 

Remark 7.4. Due to asymptotics ( 13.61 ), the condition ( 17.7b is satisfied for linearization at solitary waves. 

Proof. Without loss of generality, we will assume that ujf, = to. We have: JL{iOj)C,j = XjCj, Cj G L 2 (R n , C 2N ), and 
without loss of generality we assume that ||Cj||i 2 = L We write: 

(D m - Ui + JAj ) Cj = -V(uj (7.8) 

Let be orthogonal projections onto eigenspaces of J corresponding to ±i G c(J). Applying II to ( 17.8b and 
denoting Cf = H^Cj, we get: 

(D m - + iA,)C+ = -n+V(wj)Ci, (Pm - Uj - iAj)CJ = -irV(^)^. (7.9) 
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1A J 



Figure 1: The closest point from <r(D m — uij) to iA., is m — Uj. 



Since uij —> u>b = m, without loss of generality, we can assume that uij > rn/2 for all j G N. Since the spectrum 
cr(JL) is symmetric with respect to real and imaginary axes, we may assume, without loss of generality, that Jm\j > 
for all j € N, so that Rc iXj < (see Figure[TJ. At the same time, since Xj — > 0, we can assume that \Xj \ < m/2. 
With D m — uij being self-adjoint, one has 

|| (D m - Uij - i ) - 1 1 1 = — j- ^ rr = , (7-10) 

dist(iAj, <T{U m — uij)) \m — uij — iXj\ 

Remark 7.5. If ReiAj > to — uij, then dist(iAj, cr(D m — uij)) < \m — uij — iXj\, and the estimate ( 17.10b does not 
hold. 



Combining ( 17.91 ) and ( 17.101 ), we get 



\U-VCjW ^ C{m 2 -J]) 



Zj ls<t^ =^<^ ^r. (7.11) 

J \m — uij — iXj | \m — uij — iXj \ 



We used the normalization \\Cj\\ = 1 and the bound ||n ~V(ujj)\\L 2 ^L 2 < C\uij — u>t,\ (cf. (17.7| i). At the same time, 

*j Hz, 2 ' 



due to Re X 3 ^ 0, LemmaEUlvields = JC,-) = i||C-l| 2 - i||Cp |2 



then (17.111 ) yields 

\m — Wj — i\j\ < V2C(m 2 - uj 2 ), 

leading to 

| Xj | < V2C*(m 2 - w? ) + |m - Uj\ < (V2C + ^-) (to 2 - w?). □ 
Lemmas |6*72"1 16.31 17721 17 . 3 1 17 .41 complete the proof of Theorem |2.3l 



8 Bifurcations from the origin 

Here we prove Theorem l2.4l which we rewrite for the operator JL(w) which was defined in (14.31 ). 

Lemma 8.1. Let n < 3, and let f G C°° (R), /(ry) = r/ k + o(rj k ), fc G N. If n = % further assume that k = 1. Lef 
4>u,e~ lut be solitary waves from Lemma \3~5] Let uij — > m, ant/ Zef Aj G (T p (JL(a;j)), Aj = 0(m 2 — ui 2 ). Assume that 
the operator L_ introduced in (12.71 ) satisfies cr p (l- ) = {0}. 

7. AssMme f/iaf A = is an eigenvalue of )\ of algebraic multiplicity 2 + 2n (either n = 1, k ^ 2; or n = 2, 
k > 1; or n = 3). Then Xj =/= 0/or j G N implies that 

A b := lim — — f ^ 0. 

j'-s-oo m z — uij 

2. A b G dp Ql) or A b G (Jp(l-). 

3. If He Xj ^ 0, j G N, f/zen A^ /ias zero Krein signature: a corresponding eigenvector z satisfies 

(z,lz) = 0. 

The rest of this section is devoted to the proof of Lemma [8Tl 
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Rescaled system 



Let Cj be the eigenfunctions corresponding to Xj. Below, usually we will not write the subscripts j. The eigenfunction 
Cj corresponding to Xj € ct p (JL(cj)) satisfies 



where 



VMS = -/(CP* JPC - 2/'(C0c|>J(C3t)P<k 



(8.1) 



(8.2) 



Let e, 



cj — y "' 2 — w |' Denote Aj = Aj-/e|. Let IT* be the projections corresponding to ±i e <x(J), and let lip, 
11^4 be the projections corresponding to ±1 € <t(|3). We denote the "particle" and "antiparticle" components of C, by 
the relations 

U P C 3 (x) = Pj{ejX, ej), H A Cj(x) = eAj(ejX, e 3 ). 
We also denote II+ILp = IT|, Il+IIyi = 11^, 



U%(x) = Pf(e 3 x, ej), K%{x) = eAf(e jX , e 3 ). 



Then 



U ir^Cj(x,u;j) = ej(T)oPf)(ejX,ej), V Il A C,j(x,ojj) = e J (D Af)(e : jX,e j ). 
Let W € C 1 (M n x (0, e Q ), End(C 2JV )) be such that V(x, u) = e 2 W(ex, e). 

l 



Applying projections lip, to (18. U and dividing by e 2 in the former case and by e in the latter, we obtain the 



following system: 



m + cj 



± iA, Pf + n±W(y, e)(P i + eAj) = 0, 



-m-w±e 2 iA J )A± + en±W(y,e)(P J + 6AJ-) = 0. 



(8.3) 
(8.4) 



We took into account the relation Dofl^ = flpDo- 



We use ( 18.4b to express 



and substitute this into ( 18.3b : 



AP? 



DoPj + dgW(Pj + eAj 

m + uj =p e 2 iAj 



Aj - A+ + , 



m + uj — e 2 iAj \m + cu 



(~r^ ± iA .) ^ + n p [ w ( p . + eA .) + e 



(D W)(P J +eA J ) 



m + w =F e 2 iA 7 



(8.5) 



= 0. (8.6) 



Denoting 



and taking the sum and the difference of ( 18.61 ), we have: 



AP, 



m + oj — e 2 iAj m + uj 



P ^ l iA,Q, + npfw ( P J + e A J ) + e (DoW )^ + eA ^ 



to + uj — e 2 iAj to + uj 



iAjPj + (n+ - n p ) W(Pj + eAj) + e 



m + w =F e 2 iAj 
(D W)(P J+e A 3 ) 



= 0, 



(8.7) 



m + uj =p e 2 iAo 



= 0. (8.8) 
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The relative compactness of the eigenspaces 
Lemma 8.2. For any s € R, the multiplication operator 

z i-» wf z = wz 

is bounded from L_ s to L 2 . 

Proof. This follows from W, V y W being exponentially decaying, as functions of y = ex, uniformly as e — > 0, due 
to Lemma 1331 and due to the exponential decay of the solution u k to ( 12.6k see IIBL83I . □ 



(DqW)Z 

to + u) ^= e 2 iA 



Without loss of generality, we assume that Im < 0. Applying |Agm75 Theorem 3.2] to ( 18.61 ), we conclude 
that, for any s > 1/2, 



IIPJIIa;., < C||W 1 (P i +eA j )|| i = < CWiPj + eAj)]]^, 
for some C < oo which does not depend on j 6 N. Let us choose the normalization of C, so that 



±||2 



2 + e 



s lit*. 



= 1, 



(8.9) 



(8.10) 



Therefore, the terms in the right-hand side of ( 18.9b are bounded in L s , and this leads to 



<c, 



(8.11) 



for some C independent of j 6 N. Using (18.10b and d8.1 lb in ( 18.5b . we also get the bound for || Aj \\ H i i (as long as 

"P~l 

is bounded in H]_ 1 (W 1 , C 2 N ) . Let s > 1 ; then the embedding H}_ 1 (M™ ) C 



e > is sufficiently small). Thus, 



L 



is compact, and there is a subsequence of 



P" 
A 



, weakly in 



convergent to some 

Hg_ 1 and strongly in L 2 . Let us prove that this limit is different from zero. By Lemma l6.1l Hn+Cjl^ — ||n_C,-||x,2, 
hence 

HP+lli, +e 2 \\A+f L2 = \\r;\\h+e 2 \\Aj\\l 2 ; (8.12) 

If PJ — >• in if 1 , then, by (18.51 . A~ -> in L 2 , hence both sides of (|8.12| i converge to zero, contradicting J8. 10b . 
Recall that <Pp = nuk, where u k is a positive spherically symmetric solution to (12.6b and n 6 C w is such that 

f3n = n, ||n|| — 1; one has ITp 
The limit system 

Let W(y) = W(y, 0). Substituting lim e ^ e~ 1/{ - 2k ^(l>p{e~ l y) = <Pp(y) = nu k (y), with n e C N , \n\ = 1, we 
derive: 



n 




n 











WP = lim W(y, ey)(Pj + e 3 A 3 



<Z>p| 2fc P-2fc|<£p| 2fe - 2 (<5 P P 



-|<?p| 2fe P-2fc|<?p 



2& , 



$P 





/c 2N 



Above, ( , ) C 2w is the inner product in C 2N . Considering d8.6b in the limit e —> 0, we have 
AP± P± 



+ 



2m 2m 



2fc p ± _ 2fc?i 



21, 



,p> n± 

' C 2JV 



± iAbP* = 0, 



(8.13) 



(8.14) 



where we used ( 18. 13b . We used the equality f3n = n which implies ITp 
get the following system: 



+ 77- ul k P - 2ku k 

2m 2m k k 



21. 



Defining Q := P + — P , we 



iA b Q = 0, 



(8.15) 
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+ iA b P = 0. 



(8.16) 



n 




n 




'0" 





= ij 







in 



-^+^--^Q-2^f(M ) p) r.°' 

2m 2m fe fc \ ' /c 2N m 

We used the relation (II + — II - ) 

Proof of Lemma O (0) 

We consider the system ( 18. 7b . ( 18.81 ) (with the expression (18.5b for Aj) as a perturbation of the system ( 18.151 . ( 18.161 ). 
By IComl 11 Lemma 4.1], we have: 

^(JL(w)) D Span{Jcp, <9 w cp, d x3 cp, (a 3 - 2wa; J J)cp ; 1 < j < n], dim^(JL(w)) > 2n + 2. 



If there is an eigenvalue family (Xj)(j € N) such that Aj ^ 0, Aj 
then we would have 

dirndl) > dim^(JL(w))L <B 



>7 



0, converging to zero as uij 



rn, 



1 > 2n + 3. 



We conclude that if dim^Ql) = 2n + 2 and Aj —> 0, then Aj = for all j sufficiently large. 
Proof of Lemma 0(0 



In the case when p(y) 



, P(y)) is identically zero, ( 18.14b implies that either iAb or — iAt, belongs to 



C 2N 



<7 p (l_), with the eigenfunction being (either nonzero component of) P* e L (R",C ). Since we assumed that 
<Tp(l_) = {0}, we conclude that A b = 0. Let us note that (18.15b implies that in this case P(y) = Uk(y)M, for some 
M e C 2N . Similarly, dOBT l implies that Q(y) = u k {y)N, for some TV G C 2Ar . 

Let us now consider the case when p £ L 2 (W l , C) is not identically zero. Taking the inner product of ( 18.15b , 



(lQ6l l with 



and denoting 



,Q 



we get the equations 
which we write as 



Lp - 2fcu| fe p + iA b q = l + p + iA b q = 0, l_q + \A b p = 0, 



l_ 



A* 



(8.17) 
(8.18) 
(8.19) 



This relation implies that Af, is the eigenvalue of the linearization of the nonlinear Schrodinger equation 

itp = -— AV> - M 2k ip, V0,t)eC, iei" 

2m 

at the solitary wave u k {x)e~ luJt , ui — — ^ (°f- (1 1 -6b , ( 11.7b ). 
Proof of Lemma O(0 

If (p, g) € L 2 {W\ C) x i 2 (R", C) satisfies d8~T9l ) andp = 



, P ) , then 

c 2JV 



p = 


pn 


Q = 






\qn 




ipn 



solve d8~T5l (EJUl. It follows that P P, Q Q satisfy 







— 1_ 



1_ <£) It£2N 





p-p 

i(Q-Q). 



= A„ 



P-P 

i(Q - Q) 



(8.20) 



23 



In the case A 6 = 0, by ( 18.191 1, z 

^2N 



satisfies Iz = 0, leading to (z, \z) = 0. Then P P = u k M, Q Q = u fe AT 

for some M, AT € C 2N . As follows from ( 18.191 l, p e kcrl + = Span{<9j-u fc ; 1 < j < n}, while q e kerl_ = 
Spanju/c}. For P and Q, we have: 



Since p = 



pn 
iqn 



, P ) , M is to satisfy ( 

' C 2N \ 



Mu k , 



M 



qn 
ipn 



Nu k . 



(8.21) 



= 0. N is to satisfy N = H + M — I1~M. 



Now let us consider the case when A/, ^ 0. Since cr p (l_) = {0}, we conclude that P = P, Q = Q. Thus, 



p = 


pn 


Q = 


qn 




iqn 




ipn 



(8.22) 

Assume that the sequence Xj G (T p (JL) satisfies Rc Xj ^ 0. By Lemma loTTl the corresponding eigenvectors satisfy 

= (Cj.Kj) = (Pj + eA J; J(P, + eA,)); 

this leads in the limit e — s- to 



By (18^221 , 



The condition d8.23t leads to 



= IIP H 



= /PJP)=i||P 



+ q)n 

i(p + q)n 



+ l|2 



illP" 



(8.23) 



(P - q)n 
-i(p - q)n 



|bf + || 9 |j 2 + 2Re(p, q) |b|| 2 + \\qf - 2Re(p, g > 



2Re(p,g). 



(8.24) 



By (18.19b , the eigenvector of jl corresponding to A& is given by z = 







' 1" 




p 


>=( 








is 




-1 




}1. 




is 







using ( 18.241 ), we compute: 

= 2iHe(p,q) = 0. 



Then the relation Iz + A/,jz = leads to (z,lz) + A&(z,jz) = 0, hence (z, Iz) = 0, completing the proof of 
Lemma lSTl and thus of Theorem l2.4l 



A Appendix: Hardy-type inequalities 

Our aim in this section is to establish the generalized Hardy estimates used in the course of the study. 

A.l Hardy inequality in one dimension 

In dimension 1, Hardy's inequality reads for an integrable function /, as 



lo x ^Jo 

The result above can be has the following generalization 



±([ X \f(t)\ 2 dt)<4 f°°\f{x)\ 2 dx 



(A.l) 



Lemma A.l (Generalized Hardy inequality). Suppose that w <G C 1 ((0, oo)) is a differentiable function with positive 
derivative on (0, +oo) and / € C^ omp ((0, oo), C). Then 



w'\f\ 2 dx<4 —\f\ 2 dx. 
1 w' 
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Proof. Without loss of generality, we can assume that A = 0. The integration by parts yields 



w'(x)\f(x)\ 2 dx = - I w(x)f(x)f'(x) dx - 
Using the Cauchy-Schwarz inequality, we get: 

w'(x)\f(x)\ 2 dx<2J I w>(x)\f(x)\ 2 dx* 



w(x)f'(x)f(x) dx. 



w{x) 2 
w'(x) 



\f'(x)\ 2 dx, 



□ 



and thus / R+ w'(x)\f(x)\ 2 dx < 4 / R+ ^\f(x)\ 2 dx. 
Let U% and U\ be defined by 

U% := {ue L 2 (R,C 2 ) ; e v u e L 2 (R, C 2 )} , ^:={«e^ ; iue^}, 

endowed with the associated weighted L 2 and H 1 norms. The notation ip' stands for the radial derivative of ip. We 
will use the following elementary estimate. 

Lemma A.2. Let N = 21, I € N. Let m > anc/ fe? a, (3 be self-adjoint matrices such that a 2 = (3 2 = In, 
a/3 + j3a = 0. Then 



ii ia(m/3-A)x|| ^ r i\\ ._ / |A|+TO 

lie ii^^c^.-y^-— , 



AeR \ [-m,m] 



There is x — x{m) > such that 



ia(m(3-J)i|| 



< x(x), x € 



AeR\ (-to, to). 



(A.2) 



Proof. Let us notice that the matrix ia(m/3 — A) is skew-adjoint. Its eigenvalues are purely imaginary as long as 
A G R\[— to, to]. The norm of e 1Q ( m ^- A ) i s larger than 1 because the projections onto eigenspaces are not orthogonal. 
There are 2 x 2 Jordan blocks when A = ±to, leading to the factor of (x) in ( IA.2I ). 

Let us derive more careful estimates. Without loss of generality, we may take a = (3 = a%\ then 



M (A) := ia(m/3 - A) 



to + A 
to — A 



The eigenvalues of M(X) are ±i£(A), where £(A) = y/X 2 — m 2 > 0. For |A| > to, the operator e M W X is the 
operator of multiplication by a matrix-valued function, with the C 2 — > C 2 norm bounded uniformly in x. Without loss 

1 

of generality, let us assume that A > 0. Then the eigenvectors of M (A) are given by u± 



±i 



A — m 
A+m 



We note 



that for A 



sup || e 



with both eigenvalues of of magnitude 1, 



M(X)xn \e M ^ x {u + v)\ 

( ' Hc^t? < sup sup J : — —. — - < sup, 



\u\ 2 + \v\ 2 +2\(u,v)\ 



X u,v 



U + V 



< sup 



l(»MI 
\u\\v\ 



\u\ 2 + \v\ 2 -2\(u,v)\ - U X\ 1-M' 

\ MM 

with sup u v taken over all pairs of eigenvectors corresponding to ±i£(A). Since t^tt^t = \ — = X' one nas 

l^-HI — I 1-r rx™ 



sup 



|| e iM(A)x|| < / 



A + TO 

A — TO 



One similarly derives the estimate in the case A < — m. 

In the case when A = ±to, M(A) is a Jordan block corresponding to the eigenvalue zero; then for |A| > to and 

any x <E R one has ||e~ M ( A ) x |jc2^c 2 < with some x — x(m) > 0. □ 
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Let Q denote the operator of multiplication 

Q : u{x) i — y xu(x), u € S'(R). 

Lemma A.3. Let tp : R — > R + be a radially even differentiable function with non-negative derivative in the radial 
variable. 

IfX el\ [—to, to] and if u £ L 2 (R, C 2 ) and (D m — X)u belongs to then u G and 

\\{Q)-hi\\ ul <c{\)\\{D m -\)u\\ H o p . 
IfX G R \ (— m, m) and i/u e L 2 (R, C 2 ) aw<i (Q)(D m - X)u belongs to then u £ and 

\\(Q)- 1 u\\ n <c(X)\\(Q)(D m -X)u\\ n . 
Proof. First let us give the proof for the case when u £ C 1 (R, C 2 ) has compact support. Define 

f(x) := ia(— iud x + mj3 — X)u(x); 
one has / g (R, C 2 ) due to the compactness of its support, and we can write 

d x u + M(X)u = f (A3) 

with 

M(X) := iafim — iaX. 
The function u, being a solution to 1A.3K could be expressed as 



If | A | > to, multiplying ( IA.4I ) by the weight e v and using Lemma lAT2l leads to 
since f(x) < <p(y) for |x| < \y\. Above, Ci(A) = \ |rnr~~~ was introduced in Lemma [AT2l From the generalized 



(A.5) 



|A|-m 

Hardy inequality (cf. Lemma lATl ) with w(x) = x, it follows that 



3 2v(x) \u{x)\ 2 dx < C\(X) 2 f (f e^\f(y)\dy) dx < 4Ci(A) 2 f \x\ 2 e 2v(x) \f(x)\ 2 dx. 

JR \J-oo J JR 



By Hardy's inequality ( IA.1I ). 

I e 2v(x) \<x)\\ x < Ci[x? I l(r +a0 e <P(v)\ mdy \ 2 dx < 4Cl (Xf [ e 2 ^\f(x)\ 2 dx, 

JR+ \ x ) JR+ \ x ) K Jx ' JR+ 

and similarly for the integral over R~ . Thus, 

\\{Q)- l u\\n% <2Ci(A)||/||«o. (A.6) 

BydA3J,^((g)- 1 u) = [d x ,{Q)- 1 ]u+{Q)- 1 Mu+(Q)- 1 f- since || (Q) 2 [0x, <Q) _:l ]IU- < 1 and || A/(A)|| C 2^ C 2 = 
KWI < |A|, 

ll^((Q)-Mll«« < (|A| + l)\\(Q)-Mn% + ll<Q} -1 /llw°- 
Due to Q, \\d x ({Q)- l u)\\ n o < (2(|A| + l)d(A) + l)||/|| w o, hence 

IKQrMn < (2(|A|+2)d (A) + 1)11/11*0. (A.7) 
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If |A| > to, multiplying ( IA.4I ) by the weight and using Lemma lAT2l leads to 

\/ x ^{i-^l/CyJIdy, z>0 
Taking into account that on each of the regions of integration one has \y\ > |x — y\, one obtains : 

' ( n ~ \L + (y)\f(y)\dy, x>o 

Similarly to the previous from the Hardy inequality, there is c > such that 

WiQr'uWno <c\\(Q)f\\ K . 

Now let u £ L 2 (R, C 2 ) and (Q) (D m - X)u £ U% (E, C 2 ). We will focus on the case |A| > to; the proof in the 
case |A| > m is similar. 

First, assume that tp is bounded. If the assertion is true for C 1 functions, it is also true for square integrable function. 
Indeed, let 77 be the characteristic function of the unit ball, and set u~ £ (£) = 77(e£)u(£), e > 0; then u e £ C 1 (R) such 
that (D m - X)u e £ n°(M) and (D m - X)u £ (D m - \)u in H°(R) as e -> 0. 

In order to prove the assertion for C 1 functions, let us take a smooth cut-off function 77 which is symmetric, satisfies 
< 77 < 1, with support in \x\ < 2, and is equal to 1 for |sc| < 1. For j £ N, define r)j by r\j (x) :— r)(j). Multiplying 
by rjj the relation 

d x u + M(X)u=: f £n° v (R,C 2 ) 

gives the following relation for Vj = rjju: 

d x v 3 + M(X)vj = 77,/ + {d xrij )u £ H° V (R,C 2 ). 
Then, by (IA.7b . there is c(A) > such that 

WQ^VjuW-h^ < c(X)\\r}jf+ {d^uW-HO < c(A)(||/|| w o + i||Vi7||£~ \\X[j, 2j]u\\u%) , 

as u is in L 2 and <p is bounded, we conclude by the dominated convergence theorem in the case when cp is bounded. 
Now assume that <p is not necessarily bounded. Fix xo > 0. For j £ N, define 

<Pj(x) = if{x ) + I oX 2 ay. 



xa fa 



v>'(y) 2 +y 2 



Then ipj /• <p, ip'j — > ip', pointwise as j —> 00. Moreover, |^(a;)| is bounded, so that ||u||^o < c(A)||Q/||^o ; then 
by the Fatou lemma for the left hand side and the dominated convergence theorem for the right hand side we conclude 
that 

\\u\\ K <c(X)\\Q.f\\ n o. 

The estimates for the derivatives follows from (IA.3b. □ 



Lemma A.4. Let ip : R — > R + be a radially even differentiable function with derivative in the radial variable ip such 
that 

(v') 2 < to 2 - A 2 . 

If X £ (—to, m) and if u £ L 2 (R, C 2 ) is such that (D m — X)u belongs to H^, then u £ and 

UQr'uW^ <c(X)\\(D m -X)u\\ K . 
If X £ [—7n, to] and if u £ L 2 (R, C 2 ) is such that (Q) (D m — X)u belongs to nL, then u £ HfL and 

\\(Q)- 1 u\\ Hl <c(X)\\(Q)(D m -X)u\\ n o. 
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Proof. The proof is almost the one of Lemma |A~31 Starting from ( IA.3b . 

d x u + M(X)u = f. 
If A G (—to, to), diagonalizing the matrix M(X) gives 

dxV^ =F \/m 2 — Xu ± = f ± . 

Then 

ju-(x) = e -V m"-A"(x -,) / -( y ) ^ 

= - ; x +0 ° e v^^(x- S /) /+(y) dj/ . 

multiplying it by the weight e v provides 

[e*W\vr{x)\ < f*^ e^-^-^^^-^\e^f~{y)\ dy, 
\e^\u+(x)\ < -f+^e^-^+^^^-y^e^f+iy^dy, 

since <p(x) — \/m 2 — \ 2 x < (p(y) — \Jm 2 — \ 2 y for y < x and (p(x) + V m 2 — X 2 x < <p(y) + \Jm 2 — \ 2 y for x < y. 
The rest of the proof goes the same way as before. □ 

Remark A. 1. The above lemma can be localized outside any ball so that the assumptions can be just asserted only at 
infinity. 

A.2 Hardy inequality in higher dimensions 

Following Berthier and Georgescu |BG87|, we have the following result. 

Lemma A.5. Let n > 2. If X € C \ ±m and s > — ^ then there is a constant c — c(A, s), locally bounded in A and s, 
such that for every u G <S"(R™, C ) having the property u G Lj oc (R. n , C w ) the following inequality is true; 

\\{Q) s u\\ H i <c\\(Q) s+1 (D m -\)u\\. 

Proof. Notice that for |A| < m or A G" R, the statement is immediate due to the invertibility of D m — A. For 
A G (—00, — to) U (to, 00), the lines of the proof follow the corresponding argument by Berthier and Georgescu. Since 

\\D a (Q) s u\\ < ||[A>, (Q) s ]u\\ + ||(Q) S AHI < SlsWKQY-'uW + \\(Q) s (D m - X)u\\ + \\( m p - \)(Q) s u\\, 

it is enough to prove 

\\(Q) s u\\ <c\\(Qy +1 (D m -X)u\\. 

Using the Fourier transform, it is equivalent to proving that for every v G S' (M. n ,C N ) having the property v G 

Ll oc (R n ,C N ) one has 

\\v\\ s < c\\{h m (£) - A)u|| s+ i, 

where h m (£) is the symbol of D m . Then, up to a diagonalization, it is enough to prove that for each u G ^"(IR' 1 , C^) 
with the property u G L\ oc (W l ,C N ) one has 

\\v\\ s <c\\We + m 2 ±\)v\\ s+1 . 

Then one can either notice that the proof by Berthier and Georgescu will work, or one can recognize in the radial 
direction the symbol of the one-dimensional Dirac operator and use Lemma lA3l with ip = (s + ^i^) log(r). □ 

We can also include the thresholds at some price: 

Lemma A.6. Let n > 2. // A G C and s > — |, then there is a constant c — c(A, s), locally bounded in A and s, such 
that for every u G S"(R", C N ) having the property u G L 1 1 oc (M n , C N ) the following inequality is true: 

\\(Q) s u\\ Hl <c\\(Qy +2 (D m -\)u\\. 
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Diagonalizing the operator J, one can immediately obtain the following result: 

Lemma A.7. If A € C\ {±i(m ± cj)} and s > — i, then there is a constant c = c{X 1 uj 1 s), locally bounded in ui, A 
and s, such that for every u £ C 2N ) having the property u £ L 1 1 oc (R ra , C 2N ) the following inequality is true: 

\\(Q) s u\\ Hl < c ||(Q) s+1 (j(D m - uj) - X)u\\. 

If A £ C and s > — ~, f/ien f/;ere is a constant c — c(A, w, s), locally bounded in u), A anof s, swc/z that for every 
u £ S"(R", C 2JV ) having the property u £ L 1 1 oc (R™, C 2W ) the following inequality is true; 

\\(Q) s u\\hi < c\\(Q) s+2 {J(D m -u) - X)u\\ 

B Appendix: Carleman-Berthier-Georgescu estimates 

Lemma B.l (Carleman-Berthier-Georgescu inequality, n = 3; Theorem 5 in [BG87)). Let n = 3. For any |A| > to, 
we have the following estimate for some C, R > 0: 

Vr>l, ^||e T l a M| + ||e T l^ Ll o/3<C'||v^e^lp m -AH|, u e H^(Q R ,C N ). (B.l) 

Above, 

= {x £ K" ; |x| > i?}. 

Our aim is to have the following Carleman-Berthier-Georgescu inequalities for any dimension. 
Let D m = -iaV + /3m, 93 € C 1 ^), and denote = e v D m e- v = An + iaV^. 

Lemma B.2 (Lemma 3, IBG87I "). For v £ H^(O), 

Be((D m - iaVp + X)v, (D& - X)v) = || Vi;|| 2 + (v, [to 2 - A 2 - (V^) 2 ]i;). (B.2) 

Let r = \x\, X — xV, 3 — \{x, -iV} = -iX - f. 

Lemma B.3 (Lemma 4, BBG87I1 ). Let £ Q, <p £ C 2 {fl). Then for any u £ H^{f2), 

2 Re((£> + 2iA^ + {9, aV^})v, (Df n - X)v) 

= 2\\Vv\\ 2 +4Re(Xu, VipVv) + 2 Be(Xv, Atpv) + (v, [X(Vtp) 2 ]v). (B.3) 

Proof. First, from [$t, Do] = [— ia;V, — iaV] = [aV, uV] = aV = iZ^o, we get 

2 

4Jm{9v,(D^-X)v) = j{v, [@, D ]v) + 4Tm{@v, (Pm + iaV <p - \)v) = 2{v, D Q v) +4Re(^w, aVipv). (B.4) 
By using the identity 

(aS)(aT) = ST + i£(S,T) := ST + i^ ]k S 3 T kl S, T £ C™, (B.5) 

where the matrices = [a 3 ', a fe ] are hermitian for each j, k, we have 

D a {Df n — A) = —A + mDoP + iD (aVip) - XD = -A + mD p + Atp + iS(V, Vip) - XD Q . 

But {D Q , P} = and S(V, V<^) = S ife 9j o <p k = + ^jkfkdj = -E jk (pjd k = -S(V^, V), hence 

2 Re(A)V, (A£ - A)w) = 2|| Vw|| 2 + (w, zV) - 2i(u, £(V<^, V)v) - 2A(u, D d). (B.6) 

Adding (IB.4b and (IB.6b . we obtain 

2Re(Do«, (A™ - A)«> + 4XJm(9v, (D£ - X)v) 

= 2\\Vv\\ 2 + (v, (Aip - 2i£(Vp, V))v) + 4ARe(^f, aVtpv). (B.7) 
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Recalling that [9, Dq] = \Dq, we derive the identity 

aV<p}D = {9, aVipD } + aV^[f, D ] = {9, aVi^Do} + iaVtpD 

= ^aVrt} = -i{9 + i V^V} + {$> + ~ E(Vv?, V)}, 
where in the last line we used ( IB. 51 ). Due to J\c{S>, £(V<p, V)} = 0, the above relation leads to 
2 Re{^, qV^JDo = Im{ 2 ^ + i, V<^V} + Rc{2^ + i, S(V<^, V)} = Jm{29 + i, V^V} + 2i£(V^, V). (B.8) 
The first term in the right-hand side of ( IB. 8b is given by 



Tm{23> + i, Vv?V} = -if $>VcpV + V o Vcp9 + V o V<^ + ^VipV + iV^V - iV o Vipj 

= {Vip, V}} -Aip = 2 Tm(@{V<p, V}) - Acp = -2 Re(V o x{ Vip, V}) - Acp 

= -4 Rc(V o x VipV) - 2 Re(V o xAcp) - Aip. (B.9) 

From dOt and (|BT9T > we obtain 2 Re{^, aV^}D = -4 Re(V o a; Vt^V) — 2 Re(V o xAcp) — Acp + 2i£(V</?, V), 

2He({9,cxVcp}v,(D^- X)v) 

= 2 Re({^, aV(/)}u, £> v) + 2 Re({^, aV^}u, iaVyw) - 2A Re({^, aV^})J, «) 
= 4 Be{Xv, VipVv) + 2 Re(Xu, Aipv) - («, (Ap - 2iE(V^, V))v) 

+ (v,X(Vtp) 2 v) - A\Re{@v,aV<pv). (B.10) 

Above, we used the identity 

2Re({^, aVip}v 7 iaVtpv) = (v, ({f , a VrfiaVy - iaV(/3{^, aV^})«) = (v , X (V ip) 2 v) . 

Adding dBTTb and dBTTOt yields (|B3T >. □ 

Now we assume that ip is radial. Denote a = r~ 1 ax. We subtract ( 1B.21 > from ( IB.3b with the aid of the identity 



in 



aiTiving at 



2Re( 



{9, ot'Vcp} = {— iX — — ,acp'} = —2\acp'X — larcp" — inacp' , 



D + 2WQ) + {9, aVcp} - \(D m - iaVcp + A) v, - \)v 
( - iaV + 2(A - iacp')X + (n- |)A - f /3 - i(n - ±)<V - iarcp" v, - \)v) 



\Vv\\ 2 + A(Xv,^Xv) +2Re(Xv,Aipv) + (v 



V - m 2 + cp u + 2rcp'cp'' 



(B.ll) 



Since 



2~Re(Xv, —v) = (Xv, —v) + (v, —Xv) = (vA-Xo — -n— + —X)v) = -(v, \<p"v + (n-l)— v), 

ty rrt y> Iff an ' * f 



which is valid for v £ Hi (]?) with Q $ 0, we have 



2Be(Xv,Apv) = 2 Be(Xv, <p"v) - («, (n-l)<p" + (n - if — 



We use the above relation to rewrite ( IB .lib as 



Re 



- iaV + 4(A - \atp')X + (2n - 1)(A - iaip') -m/3- 2\arip" v, - X)v 



II Vw|| 2 + A\\(^)iXv\\ +2Be(X Vl y"v) + (v, A - m 2 + ip' 2 + 2r V 'y" - (n - 1)^" - (n - l) 2 ^ 
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For any S > and any j& (r) > 0, the above relation yields the following inequality: 



[ - iaV + 4(A - ]&<p')X + (2n - 1)(A - ia<p') -mp- 2ianp"] - — .„ 
> \\Vv\\ 2 + {A-8)\\{^-) x l*Xvf - S^Wi^-fiyvf 



+ \\^ 2 (D* l -\)vf 



+(v, A 2 - to 2 + ip' 2 + 2r(p'ip" - (n - l)<p" -(n - l) 2 y v}. 
To eliminate from (IB. 12b the terms with Vu and Xv, we need to take 5 > and ^(r) so that 
i— [ - iaV + 4(A - xay')X]u\\* < II V«|| 2 + (4 - SM^y^XvW 



(B.12) 



"2^/2 
which would follow from 

It is enough to take 

Now JB. 12b takes the form 
1 



— <1, 


A 2 + v ' 2 , 


r 


<S = 2, 


#(r) = ^ + 2^I- 
4 93' 





(B.13) 



|| . (2n - 1)A - to/9 - i(2n - 1)<V - 2iar<p" t>|| 2 + ||^ 1/2 (X>£ - A)w|j' 



2.^V2 

> (v, A 2 - m 2 + y/ 2 + 2rvV" - (n - 1 V - (n - 1) 

which we rewrite as 

W^t^D* - A) V || 2 > /«, [A 2 - to 2 + <^ 2 + 2rvV 



2^ ^ 



//2 



(B.14) 



-(n- %" - (n-1) 



2 93' V 2 (2n-l)(|A|+¥> / )+TO + 2r|v3 // | 



«). (B.15) 



r 293' 4^#(r) 

Lemma B.4 (Carleman-Berthier-Georgescu inequality in R"). Lef A G R. Let if € C 2 (R + ) vWf/z > 0. Asswme 
f/iflf f/;ere is i?^ > such that for all r > one has 

J2 



A 2 - m 2 + %- > 0, 



1/9 „ / ;/ / ^ // / rco" 2 3n to |w"| 

+ 2W > (» - + (n - lf V - + + ^ + ^ + ^ 



(B.16) 
(B.17) 



Then 



X 2 - m 2 + ^-Y"e^v\\ < \\(~+2^-+2r<p'y / *e'0(D m -\)v\\, v e Hl{Q Rv ) 



1/2 



,n A r 



V2, 



"( 



-^ + 2V) 1/2 ' 
-' 2 \V2 ,„ „ _ „ ,„ A 2 r 



Moreover for any v £ Hl(Q Rtf ) such that (f + 2^f + 2r</) ' e v (D m - X)v £ L 



2/inm <rN\ 



A 2 - to 2 + fO\ V >„|| < II (~ + + 2V) V Vp m - A)«| 



Proof. One has 



(2"-l)|A| < (2n-l)|A| 



n+gi^+SrV — 16r|A| 8r<p' — 8r 



< hence the inequality ( IB. 17b implies that, as long as r > i?^, 



the following relation holds: 

V + 2W > (» - + (n - l) 2 ^ + ^ + ^ - 1)(|A| + ^ + " + 2 ^"' 
2 r 2</?' 



n + + 8r<^ 



(B.18) 



Now we use the inequality (IB. 15b (with given by ( IB. 13b ) together with (IB. 18b . substitute e v v in place of v, and 
use the identity D l ^ l (e v v) = e v D m v, arriving at the desired inequality. 

The extension to non compactly supported functions follows exactly the process of the proof of Lemma lA31 □ 
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